.  J  :T.‘  i  1 


.si* -‘A  -  ''iy> 

■  *»  "  *  •  r 

I'  ;  .. 

ARLTR-86-14 

■  '  .  ' 

!  i ■  "  r 1  ' 

*  .  J 

Copy 

AN  EXPERIMENTAL  INVESTIGATION  OF  ACOUSTIC  PROPAGATION 
IN  SATURATED  SANDS  WITH  VARIABLE  FLUI&  MMttfmtS 


O 

00 

R.  Daniel  Conley,  Jr. 

00 

cr> 

LO 

« — 

APPLIED  RESEARCH  LABORATORIES 

< 

i 

Q 

THf  UNIVERSITY  OF  TEXAS  AT  AUSTIN 

FOOT  OFFICE  SOX  SOM.  AUSTIN,  TEXAS  TBTIT  S01S 

< 

May  1985 

Technical  Report 

Approved  for  public  ralaaaa; 

dittributlon  unlimited. 

Prepared  for: 


G"3 


OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VA  22217 


DTIC 

elect* 


W 


OCT  9  1905 


A 


85  AO  9  083 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

!.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO 

4M  /.■>■“?/ 

3.  RECIPIENT'S  CATALOG  NUMBER 

■id 

4.  TITLE  (mid  Subtlllo) 

AN  EXPERIMENTAL  INVESTIGATION  OF  ACOUSTIC 

5.  TYPE  OF  REPORT  6  PERIOD  COVERED 

PROPAGATION  IN  SATURATED  SANDS  WITH  VARIABLE 

FLUID  PROPERTIES 

6.  PERFORMING  ORG.  REPORT  NUMBER 

ARL-TR-85-14 

7.  AUTM0R(»J 

R.  Daniel  Costley,  Jr. 

8.  CONTRACT  OR  GRANT  NUMBERf*) 

N00014-80-C-0490 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Applied  Research  Laboratories 

The  University  of  Texas  at  Austin 

Austin,  Texas  78713-8029 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 

AREA  a  WORK  UNIT  NUMBERS 

technical  report 

II.  CONTROLLING  OFFICE  NAME  AND  ADORESS 

Office  of  Naval  Research 

12.  REPORT  DATE 

May  1985 

Department  of  the  Navy 

Arlington,  Virqinia  22217 

14.  MONITORING  AGENCY  NAME  8  AGDRESSf//  dllltrent  Iron  Controlling  Olllco) 

13.  NUMBER  OF  PAGES 

100 

15.  SECURITY  CLASS,  (o(  thlt  ropnrt) 

UNCLASSIFIED 

1S«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  thlt  Roporl) 

Approved  for  public  release;  distribution  unlimited. 


18.  SUPPLEMENTARY  notes 

.....  A  Kj 


ELECTE] 
OCT  9  1985 


19.  KEY  WOROS  (Continue  on  reverae  aide  It  neceaaery  and  identity  by  block  number) 

^ompresslonal  velocity-  experiment 

compress ional  attenuation-  sand 

ITiot  theory; 
viscosity  dependence  - 

20.  ABSTRACT  fConfinu*  on  reverae  aide  If  neceeaery  end  Identity  by  block  number) 

?The  Biot-Stoll  theory  describes  the  propagation  of  acoustic  waves  in  a  saturated, 
unconsolidated  porous  medium.  The  expressions  for  the  attenuation  and  phase 
velocity  derived  from  this  theory  depend  explicitly  on  the  viscosity,  density, 
and  bulk  modulus  of  the  pore  fluid.  An  experiment  has  been  designed  to  determine 
[  the  dependence  of  attenuation  and  phase  velocity  on  these  properties  of  the  pore 
fluid.  The  phase  velocity  and  attenuation  of  compressional  waves  were  measured 
using  a  mixture  of  water  and  glycerine  as  the  Interstitial  fluid.  The  theore¬ 
tical  background  is  reviewed  and  the  experimental  procedure  is  discussed  in  — ^ 


FORM 
1  JAN  73 


EDITION  OF  1  NOV  65  IS  OBSOLETE 


UNCLASSIFIED^ 

SE CURITY  CLASSIFICATION  OF  THIS  PACE  'When  Data  Enfrtdf 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGCfWfe«t  Dmtm  Entmfd) 


20.  ^(cont'd) 

tail.  The  results,  along  with  comparisons  with  the  Biot-Stoll  theory, 
are  then  presented.  The  choices  of  the  theoretical  parameters  are  discusse 
and  their  relation  to  the  fit  of  the  theory  to  the  data.  The  Biot-Stoll 
theory  is  shown  to  adequately  describe  the  effects  of  the  fluid  properties 
on  acoustic  wave  propagation  in  saturated  sediments,  at  Ijeast  for 
compressional  waves  of  the  first  type.  . /r'^yc^.ws  •  ,  ■  /•' 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  P  AG  Ef  *»>•»*  Enffd) 


'pv 


TABLE  OF  CONTENTS 


Page 


LIST  OF  TABLES  v 

LIST  OF  FIGURES  vil 

I.  INTRODUCTION  1 

II.  THEORETICAL  BACKGROUND  4 

III.  LABORATORY  EQUIPMENT  11 

A.  General  Procedure  1 1 

B.  Instrumentation  13 

C.  Transducers  15 

1.  Construction  15 

2.  Geometrical  Spreading  17 

3.  Directivity  20 

D.  Sediment  Preparation  24 

IV.  EXPERIMENTAL  PROCEDURE  29 


A.  Phase  Velocity 

B.  Attenuation 

1.  Measurement  Technique 

2.  Comparison  of  Results 

3.  Transient  Consideration^ 

4.  Effects  of  Dispersion 


Accesion  For 

NT  1 3  CRA&I 
OTIC  TAB 
U:  cni.ced 

J  .li  c.it, on 


Dy 

D. 


Di.t 


;  or 

t>o  .  jl 


29 


35 


35 

44 

47 


54 


Page 


V.  EXPERIMENTAL  RESULT  56 

A.  Physical  Properties  of  the  Sediment  56 

B.  Comparison  of  Results  62 

VI.  DISCUSSION  AND  CONCLUSIONS  69 

APPENDIX  A  -  COMPUTER  PROGRAM  72 

APPENDIX  B-  EXPERIMENTAL  DATA  79 

REFERENCES  95 


w 


5^®2S^SPS9BSC^BS058H®SBWS5RPJ'3SKB!^C!KE*3!^SnWOT!3!SC3BC3*3Jr^CTOI^t  ^/n?ntr^ia~*,'c*M'j^i  xrK*r%:Kr^j*x-w*'i™  xru  i."*  a  -u-.^n*A"«!«T(-irw4i-"*t 


LIST  OF  TABLES 


Table  Page 

I  Comparison  of  Velocity  Results  33 

in  MS-MH  and  MS-L  Size  Beads 

II  Analysis  of  Errors  in  Attenuation  41 

III  Comparison  of  Attenuation  Results  46 

in  MS-MH  Size  Beads 


IV  Comparison  of  Attenuations  and  Standard  Deviations  55 
Measured  from  Different  Cycles 
of  the  Same  Received  Signal 


V  Physical  Properties  and  Parameters  of  Sediment  57 

Used  in  Theoretical  Calculations 

VI  Physical  Properties  of  Fluid  Used  59 

in  Theoretical  Calculations 

Vh  Velocity  and  Attenuation  in  MS-MH  Beads  64 

with  Variable  Fluid  Properties 

VIII  Porame.ers  Used  in  FORTRAN  Program  6IV0  73 


v 


I.  INTRODUCTION 


There  has  been  continued  interest  in  acoustic  propagation  in 
fluid  saturated,  porous  media  for  the  past  three  decades.  Initially,  much 
of  the  interest  was  due  to  problems  encountered  in  the  oil  industry. 
Acoustic  techniques  were  devised,  and  are  still  being  developed  and  used, 
to  determine  the  properties  of  the  porous  rocks,  and  their  interstitial 
fluids,  found  at  different  depths  in  oil  wells.  Maurice  Biot,  while 
working  for  Shell  Development  Company,  published  two  papers  in  1956 
end  a  third  in  1962  describing  wave  propagation  in  fluid  saturated, 
porous,  elastic  solids.  This  theory  was  used  to  model  wove  propagation 
in  fluid  saturated,  porous  rocks  where  the  individual  grains  are  cemented 
together. 

Since  then,  there  has  been  interest  in  accurately  describing 
wave  propagation  in  unconsolidated,  fluid  saturated,  porous  medio,  such 
os  water  saturated  sands.  This  topic  has  widespread  application  in 
marine  seismic  surveys  and  in  sonar  and  underwater  acoustics,  since 
portions  of  the  ocean  bottom  consist  of  water  saturated  sand.  It  is  also 
of  interest  to  soil  mechanicists,  who  use  seismic  methods  to  determine 
the  structural  properties  of  sends  and  soils,  both  onshore  and  offshore. 

To  predict  the  acoustic  properties  of  these  unconsolidated 
sediments,  R.  D.  Stoll  extended  Biot’s  original  theory  to  take  into  account 
the  losses  at  the  grein-to-grein  contacts  by  proposing  that  some  of  the 
constitutive  coefficients  in  Biot’s  equations  be  made  complex  (Stoll, 
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1974).  Thus,  he  modeled  the  solid  matrix  as  a  viscoelastic,  porous  solid, 
rather  than  as  the  porous,  elastic  solid  considered  in  Biot's  original 
work. 

Due  to  the  difficulties  encountered  when  making  measurements 
in  saturated  sediments,  only  a  limited  amount  of  date  exists  in 
publication.  Some  published  data  on  unconsolidated,  fluid  saturated, 
porous  media  have  been  produced  by  R.  D.  Stoll  (Stoll,  1974),  David  Bell 
(Bell,  1979),  D.  J.  Shirley  (Shirley,  1978;  Bedford,  1962),  end  Jens  Hovem 
and  G.  D.  Ingram  (Hovem,  1979).  More  measurements  need  to  be  token, 
however,  both  to  validate  the  theory  and  to  determine  the  physical 
properties  of  the  sediments. 

The  present  experiment,  first  proposed  by  D.  J.  Shirley, 
investigates  the  dependence  of  the  attenuation  end  phase  velocity  of 
compressionol  waves  in  a  saturated  sand  on  the  properties  of  the 
interstitial  fluid.  A  mixture  of  glycerine  end  water  was  used  os  the 
saturating  fluid  to  permit  the  controlled  variation  of  the  viscosity, 
density,  and  bulk  modulus  of  the  pore  fluid.  Since  these  parameters 
appear  explicitly  in  the  Biot-Stoll  equations,  this  investigation  provides 
a  critical  test  of  the  theory.  Furthermore,  the  measurements  also 
provide  an  estimate  of  some  of  the  constitutive  coefficients  in  the 
equations. 

Similar  measurements  have  been  made  by  Shirley,  and  later  by 
Elliot  (Bedford,  1982).  Due  to  large  experimental  uncertainties  in  their 
results,  no  definite  conclusions  could  be  made.  After  investigating  their 
procedure,  a  few  simple  modifications  have  been  made.  New 


measurements  hove  been  token  and  shown  to  be  repeatable  within 
experimental  uncertainties. 

The  most  significant  development  over  the  earlier  work  resulted 
from  the  redesign  of  the  compressional  wave  transducers.  With  this 
improvement,  the  existence  of  spherical  spreading  was  established  at 
the  separations  at  which  these  measurements  were  taken  and  the  error 
caused  by  the  directivity  of  the  transducers  was  reduced. 

Both  Shirley  and  Elliot  made  shear  wove  measurements  in  their 
investigations.  Although  it  would  have  been  desirable  to  hove  measured 
the  attenuation  and  phase  velocity  of  shear  waves,  it  was  decided  that 
we  could  not  substantially  improve  their  earlier  results  without 
considerable  additional  investment  in  time  and  effort.  Therefore,  only 
compressional  wave  propagation  was  considered  in  this  investigation. 

This  report  is  organized  as  follows.  A  short  review  of  the 
theory  for  the  compressional  wave  case  is  presented  in  Chapter  II.  In 
Chapters  III  arid  IV,  the  experimental  procedure  is  discussed  and  some  of 
the  assumptions  are  justified.  We  compare  our  results  to  those  of  other 
investigators  in  Chapters  IV  end  V,  end  in  Chapter  V  our  results  ore 
compared  with  the  theory. 


II.  THEORETICAL  BACKGROUND 


In  three  articles  published  in  The  Journal  of  the  Acoustical 
Society  of  America,  two  in  1956  arid  the  third  in  1962,  M.  A.  Biot 
presented  a  theory  describing  the  propagation  of  acoustic  waves  in  a 
fluid  saturated,  porous,  elastic  solid  (Biot,  1956a,  1956b,  1962)  In  this 
theory,  dissipation  due  to  scatter  is  neglected,  requiring  that  the 
wavelength  be  large  in  relation  to  the  pore  size. 

Biot  predicted  three  types  of  body  waves  in  an  unbounded, 
porous  medium:  two  types  of  compressional  waves  and  one  shear  wave. 
All  three  of  these  have  been  observed  experimentally.  The  first  type  of 
compressional  wave  and  the  shear  wave  are  analogous  to  the  waves 
predicted  by  ordinary  elasticity.  The  second  type  of  compressional  wave 
is  more  dispersive  and  highly  attenuated  (Stoll,  1974,  1977). 

In  Biot's  theory,  the  solid  matrix  is  an  elastic  frame  in  which 
there  is  no  relative  motion  between  the  groins.  Later,  Stoll  and  Bryan 
extended  this  theory  to  include  unconsolidated  porous  media,  such  os 
sands  and  sediments.  Thus,  they  incorporated  into  the  theory  the 
inelastic  nature  of  the  frame,  which  is  caused  by  the  relative  motion  of 
the  grains. 

In  the  Biot-Stoll  theory,  losses  ore  attributed  to  three 
mechanisms.  There  ore  frictional  losses  at  the  groin-to-groin  contacts. 
Secondly,  the  movement  of  grain  particles  relative  to  one  another 
produces  fluid  motion  in  their  immediate  vicinity  that  does  not 
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contribute  to  the  propagation  of  the  wove.  This  local  fluid  motion  is  at 
the  expense  of  the  propagating  wave,  and  thus  it  constitutes  a  loss 
mechanism. 

The  third  type  of  loss  mechanism,  and  the  one  which  is  thought 
to  dominate  in  this  study,  is  due  to  viscous  losses  in  the  fluid  This  type 
of  absorption  becomes  important  if  there  is  significant  motion  between 
the  fluid  and  the  skeletal  frame.  This  effect  dominates  at  higher 
frequencies  or  in  sediments  where  the  permeability  is  high.  For 
instance,  this  effect  will  become  more  important  at  lower  frequencies 
for  coarse  sands  than  it  would  for  finer  sediments  and  clays,  which  have 
a  lower  permeability. 

The  main  purpose  of  this  report  is  to  compere  the  experimental 
measurements  of  attenuation  and  phase  velocity  with  the  predictions  of 
the  theory.  Whet  follows  is  an  outline  of  the  derivation  of  an  expression 
for  the  wave  number  from  Biot's  equations.  It  will  serve  only  to  acquaint 
those  already  familiar  with  the  theory  with  our  method  of  calculating 
the  phase  velocity  and  attenuation,  should  there  arise  any  discrepancies 
between  our  results  and  those  of  others.  For  a  more  complete 
discussion,  the  reader  should  consult  the  references  by  Biot,  Stoll, 
Bedford,  and  Hovem  and  ingram 

In  order  to  present  Biot's  equations,  some  of  the  notation  used 
needs  to  be  introduced  If  we  let  u  and  U  represent  the  displacement 
fields  of  the  solid  skeletal  frame  and  the  fluid,  respectively,  end  if  $  is 
the  porosity,  then  we  can  write 
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e  =  div  u_=  V  •  u 

and  £  =  £div(u_-y)  =  £V-(u-IJ) 

Thus,  e  represents  the  dilation  of  an  element  attached  to  the  frame,  and 
$  represents  "the  volume  of  fluid  that  has  flowed  in  or  out  of  an  element 
of  volume  attached  to  the  frame"  (Stoll,  1974).  The  mass  density  of  the 
fluid-solid  mixture  is 


and 


P  =  <  1  -  £ )  Ps  *  fi  Pi 
Pc  =  p,/f  *  C/y 


where  ps  and  pf  represent  the  density  of  the  solid  and  fluid. 


respectively. 

The  dreg  and  virtual  mass  coefficients  are  b  and  c,  respectively. 
Hovem  and  Ingram  (Hovem,  1979)  evaluated  b  and  c  as  funtions  of 
frequency  and  showed  that 

by  .  Tt  Real  F(x), 

cy  a  imaginary  FOc)^ 

where  k2  =  ap2  ( wPf/T )  , 

t\  represents  the  viscosity,  u>  is  the  angular  frequency,  and  B(l  is  the 
permeability.  The  term  ap,  which  represents  a  pore  size  parameter,  is 
given  by 

ap  = e  dm/  j<  i  - 1 ) 

and  is  the  diameter  of  the  spherical  grains  of  the  sediment  The 
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function  F(x)  can  be  evaluated  in  terms  of  the  Kelvin  function  T(x) 


kT(k)/4 


F(x)  = 


1  _  2  TOO/ 


IX 


The  Kelvin  function  is  given  by 


ber‘(x)  +  i  bei'(x) 


T(k)  = 


ber(x)  +  i  bei(x) 


The  functions  ber(K)  and  bei(x)  are  the  real  and  imaginary  parts  of  the 
zero  order  Bessel  function  of  the  first  kind 

ber(x)  +  i  bei(x)  =  J0  (  x- e’3TT/4  ) 

The  functions  ber'(x)  and  bei'(x)  are  the  derivatives  of  ber(x)  and  bel(x), 
respectively. 

Biot’s  equations  contain  several  constitutive  coefficients  which 
can  be  evaluated  in  terms  of  the  moduli  of  the  solid  matrix  and  the  solid 
and  fluid  constituents.  The  bulk  modulus  of  the  drained,  solid  frame,  Kb, 

and  the  shear  modulus  of  the  drained,  solid  frame,  u,  are  made  complex 
in  order  to  account  for  the  inelastic  response  of  the  frame.  The  bulk 
modulus  of  the  solid  and  fluid  constituents  are  denoted  K$  and  Kf  , 

respectively.  In  terms  of  these  parameters,  the  constitutive 
coefficients  in  Biot's  equations  can  be  written  (Stoll,  1974) 
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of  the  wove  number  by  the  following  relotions: 

Re(k)  =  w/Co  and  !m(k)  =  -« 

Substituting  equations  (11-1)  into  Siot’s  equations  leads  to  two 
algebraic  equations  in  c  and  £  which  can  be  written 

(  -o2p  ♦  k2H  )  e  +  (  o>2pf  -  k2C )  £  =  0 

(  -<»>2pf  +  k2C  )  e  +  (  -<i>2pc  -  k2M  -  )  £  =  0 

In  order  that  nontrivial  solutions  to  this  system  of  equations  exist,  it  is 
required  that  the  determinant  of  the  coefficients  vanish.  Solving  this 
determinant  leads  to  a  fourth  order  algebraic  equation  in  k 

A,  k4  +  A2  k2  ♦  A3  =  0 
where  A,  =  C2  -  hH  , 

Aj  =  (pM  »  PCH  -  2pfC  )  - '“HD/#2  _  („-2) 

and  A3  =  ( 1pu>3  )b/^2  +  u>4  ( pf2  -  ppc  ) 

The  roots  to  this  equation  are  given  by  the  quadratic  formula 

k2  =  (  -A2  t  /(A22  -  4A,A3)  1  /  2A,  .  (II-3) 

This  represents  four  roots,  two  positive  and  two  negative. 

The  two  positive  roots  correspond  to  two  different 
compressional  waves.  The  value  of  k  calculated  from  the  negative  value 
of  the  discriminant,  -/(A22  -  4A,A3),  corresponds  to  the  compressional 
wave  of  the  second  type,  or  the  compressional  slow  wave.  Since  the 
present  experiments  are  designed  to  measure  the  velocity  end 
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attenuation  of  the  compressional  wave  of  the  first  type,  the  wave 
number  was  calculated  with  the  positive  discriminant. 

The  FORTRAN  program  shown  in  Appendix  A  was  used  in  order  to 
solve  the  quadratic  formula  (11-3)  with  coefficients  (11-2)  as  a  function 
of  viscosity.  The  program  displays  the  results  in  graphical  form,  in 
order  to  solve  for  the  wave  number,  and  thus  the  phase  velocity  and 
attenuation,  the  physical  parameters  of  the  saturated  sediment  need  to 
be  determined.  This  process  is  discussed  more  fully  in  Chapter  V.  In  the 
next  two  chapters  the  procedure  used  to  measure  these  quantities  is 
discussed  in  detail. 


III.  LABORATORY  EQUIPMENT 


A.  General  Proj&dio 

The  objective  of  this  experiment  is  to  determine  the  phase 
velocity  and  attenuation  of  harmonic,  compressional  waves  in  a  fluid 
saturated  sediment  as  a  function  of  the  properties  of  the  pore  fluid. 
These  results  will  be  compared  with  those  predicted  by  the  Biot-Stoll 
theory.  Since  the  fluid  viscosity,  bulk  modulus,  and  density  appear  in  the 
expression  for  the  wave  number,  these  measurements  are  expected  to  be 
a  critical  test  of  the  theory. 

Similar  measurements  of  the  acoustic  properties  of  sediments 
have  been  performed  by  Bell  (Bell,  1979),  Shirley,  and  Elliot 
(Bedford,  1982).  The  main  focus  of  the  work  described  in  this  chapter  and 
the  next  has  been  to  improve  the  techniques  developed  by  these 
investigators  for  measuring  the  attenuation  and  phase  velocity  of 
compressional  waves  and  to  locate  the  causes  of  experimental  error.  The 
most  significant  improvement  over  the  earlier  work  has  resulted  from 
the  redesign  of  the  compressional  wave  transducers.  With  this 
improvement,  spherical  spreading  was  established  at  the  separations  at 
which  these  measurements  were  taken  and  the  error  caused  by  the 
directivity  of  the  transducers  was  reduced. 

The  two  transducers  ore  essentially  identical.  One  serves  as  a 
transmitter  and  the  other  os  a  receiver.  The  transmitter  converts  on 
electrical  signal,  which  is  generated  in  on  oscillator,  into  on  acoustic 
signal,  which  propagates  through  the  sediment.  The  receiver  detects  this 
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acoustic  signal  and  reconverts  it  into  an  electrical  signal  which  is 
displayed  on  on  oscilloscope. 

The  central  idea  in  our  measurement  procedure  is  to  compare 
the  time  delay  and  amplitude  of  the  output  voltage  signal  at  one 
transducer  separation  with  the  same  measurements  at  other  separations. 
From  the  time  delay  measurements  at  different  separations,  the  phase 
velocity  is  determined.  Similarly,  by  compering  the  different  amplitude 
measurements,  the  attenuation  can  be  calculated.  It  is  assumed  that  the 
coupling  between  the  transducer  and  sediment  is  the  some  at  each 
separation.  Thus  we  assume  that  the  decay  in  amplitude  is  due  only  to 
the  increase  in  transducer  separation. 

Initially,  the  transducers  are  placed  20.0  cm  apart.  The  time 
delay  from  the  beginning  of  the  triggered  pulse  to  the  fourth  positive 
peak  of  the  received  signal  is  then  recorded  along  with  the  peak-to-peak 
amplitude  of  the  fourth  cycle  of  the  received  signal.  The  fourth  cycle  is 
used  so  that  the  measurements  ore  made  as  close  to  the  steady  state 
port  of  the  signal  os  possible.  The  signal  is  then  monitored  for  at  least 
four  hours  to  ensure  that  there  are  no  variations  in  the  signal  amplitude 
due  to  changes  in  the  coupling  between  the  sediment  and  the  transducers 
or  due  to  any  settling  of  the  sediment.  Four  to  six  readings  are  taken 
over  this  period,  and  the  slight  variations  about  the  stabilized  value  ore 
averaged.  The  receiver  is  then  moved  2.0  cm  closer  to  the  transmitter 
and  the  amplitude  and  time  delay  measurements  ore  repeated.  The  entire 
procedure  is  repeated  at  2.0  cm  intervals  until  the  separation  between 
the  transducers  has  been  reduced  to  6.0  cm.  Using  this  set  of 


measurements,  the  phase  velocity  end  attenuation  ore  determined  by 
finding  the  least  squares  fit  to  the  data. 

The  fluid  properties  are  then  altered  by  adding  glycerine  to  the 
mixture.  The  entire  set  of  measurements  is  repeated  and  the  phase 
velocity  and  attenuation  ore  determined  for  a  new  value  of  viscosity. 
This  procedure  is  repeated  several  times,  each  time  with  a  new  value  of 
viscosity.  In  this  way,  plots  of  phase  velocity  and  attenuation  versus 
viscosity  con  be  mode  and  then  compared  with  predictions  of  the  theory. 

B.  Instrumentation 

A  schematic  of  the  laboratory  equipment  is  shown  in  Fig.  1.  A 
continuous  tone  of  105.0  kHz  and  10.0  v  peek-to-peak  amplitude  is 
generated  in  the  oscillator.  The  output  is  sent  to  a  frequency  counter  and 
a  tone  burst  generator.  The  tone  burst  generator  passes  a  four  cycle 
pulse  of  105.0  kHz  to  a  transmitting  transducer  and  an  oscilloscope.  The 
signal  is  then  cut  off  for  several  milliseconds  so  that  the  reflections  of 
the  signal  off  of  the  sides  of  the  container  and  the  sediment-water 
interface  will  die  down  before  a  new  four  cycle  pulse  is  sent  to  the 
transmitter.  If  the  reflections  are  not  suppressed,  they  will  interfere 
with  the  signal  that  follows  the  direct  path  between  the  transducers  and 
produce  an  error  in  the  signal  amplitude  and  time  delay  measurements. 

The  pulse  is  also  displayed  ori  an  oscilloscope.  The  tone  burst 
generator  has  a  synchronized  output,  which  is  connected  to  the  external 


trigger  input  of  the  oscilloscope.  This  causes  the  oscilloscope  trace  to 
trigger  at  the  beginning  of  the  four  cycle  pulse. 

The  acoustic  pulse  generated  in  the  sediment  is  then  detected  by 
the  receiving  transducer.  The  transducer  translates  the  acoustic  signal 
into  a  voltage  which  is  sent  through  a  bandpass  filter  and  then  displayed 
on  the  oscilloscope.  The  low  and  high  cutoff  frequencies  of  the  bandpass 
filter  are  100  kHz  and  110  kHz,  respectively.  The  filter  is  used  to 
remove  the  extraneous  frequencies,  introduced  because  the  transmitted 
pulse  is  the  product  of  a  pure  tone  and  the  difference  of  two  step 
functions,  from  the  received  signal. 

C.  Transducers 

1.  Construction.  The  two  transducers  are  nearly  identical  and 
can  be  used  as  either  transmitters  or  receivers.  The  transducers  were 
built  from  a  piezoelectric  element  mounted  in  a  stainless  steel  plate. 
A  rod  was  welded  onto  the  plate  so  that  the  plate  and  transducer 
element  could  be  inserted  into  the  sediment-water  mixture.  A  sketch 
of  a  transducer  mount  is  shown  in  Fig.  2. 

The  two  elements  are  made  of  the  materiel  Chennelite  5500  and 
manufactured  by  Channel  industries,  Inc.  Their  thickness  and  diameter 
are  3/16  in.  and  3/4  in.,  respectively.  They  each  have  two  resonances. 
In  the  thickness  mode,  #1  and  *2  resonate  at  446.4  kHz  and  445.2  kHz, 
respectively.  In  the  radial  mode,  #t  and  *2  resonate  at  105.5  kHz  and 
105.1  kHz,  respectively.  In  this  experiment,  the  transducers  are  used 
in  the  redial  mode,  at  105.0  kHz. 
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To  build  the  transducers,  a  7/8  in.  hole  was  drilled  through  the 
steel  plate.  The  element  was  placed  in  the  hole  with  corprene  placed 
between  it  and  the  steel  plate,  around  the  circumference  of  the  element, 
and  behind  it  on  the  backside  of  the  transducer.  A  wire  was  soldered  on 
each  face  of  the  element,  and  the  two  wires  were  placed  in  a  shielded, 
waterproofed  sheath. 

The  corprene  is  of  much  lower  acoustical  impedance  than  either 
the  piezoelectric  element  or  the  sediment-water  mixture.  It  was  thus 
used  on  the  backside  of  the  element  to  act  as  a  pressure  release  material 
so  that  most  of  the  radiation  would  be  from  the  front  side  of  the 
transducer.  This  directs  most  of  the  energy  toward  the  receiver, 
producing  a  higher  signal  to  noise  ratio. 

The  area  of  the  plate  was  made  sufficiently  large  and  thick,  arid 
the  shape  of  the  element  was  chosen  to  be  circular,  so  that  radiation 
from  the  transducer  could  be  treated  as  radiation  from  a  piston  in  a  rigid 
baffle.  In  this  way  the  geometrical  spreading  and  the  directivity  of  the 
transducer  could  be  calculated  and  verified  by  measurement. 

2.  Geometrical  Spreading.  Since  the  geometrical  spreading  in 
the  neerfield  is  very  complicated  and  varies  with  distance  from  the 
source,  the  transducers  were  designed  for  use  in  the  farfield.  The 
Rayleigh  distance  is  roughly  the  distance  away  from  the  source  at  which 
the  farfield  begins.  At  separations  larger  than  the  Rayleigh  distance,  the 
amplitude  of  the  signal  decays  as  1/r,  due  to  the  geometrical  spreading 
of  the  signal. 
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The  Rayleigh  distance  for  a  circular  piston  in  an  infinite,  rigid 
baffle  is  given  by  (Blackstock) 

R  =  V2  ka2 

where  a  is  the  radius  of  the  piston,  k=(o/c0  is  the  wave  number,  and  <.o  is 
the  angular  frequency.  Note  that  R  is  inversely  proportional  to  the  phase 
velocity  c0.  For  these  transducers,  and  a  frequency  of  105.0  kHz,  this 
formula  gives  a  Rayleigh  distance  of  2.0  cm  in  water.  Thus,  spherical 
spreading  is  expected  at  separations  of  about  4.0  cm  or  greater. 

For  a  lossless  medium  in  which  there  is  l/r  spreading,  the 
signal  amplitude  satisfies  the  following  relation: 


E  KA 


where  E  is  the  signal  amplitude,  Ef.  is  a  reference  amplitude,  A  is  the 

wavelength,  and  K  is  a  calibration  constant.  Taking  the  logarithm:  of  both 
sides  gives 

E  r 

log  —  =  -  log  —  ♦  log  K . 

Er  A 

Thus,  if  the  signal  amplitude  is  measured  as  a  function  of 
separation  for  separations  greater  than  the  Rayleigh  distance,  the  points 
will  fall  on  a  sti eight  line  with  a  slope  of  -1  when  plotted  on  log-log 
graph  paper.  Thus,  to  experimentally  verify  that  there  is  spherical 
spreading  in  water  at  separations  greater  than  4  0  cm,  the  amplitude  of 


the  acoustic  signal  was  measured  as  the  separation  was  varied  from 
4.0  cm  to  20.0  cm,  in  2.0  cm  increments.  A  least  squares  fit  through  the 
points  on  a  log-log  plot  of  range  versus  amplitude  revealed  a  slope  of 
mr- 99 1 4  of  the  equation. 

E  r 

log  —  =  m  x  log  —  +  log  K 
Er  X 

!f  the  point  corresponding  to  r  =  4.0  cm  is  omitted,  the  slope  is 
m=- 1.0061  as  shown  in  Fig.  3.  At  these  distances,  it  can  be  assumed  that 
there  is  no  absorption  in  water  and  the  total  losses  can  be  attributed  to 
geometrical  spreading.  Therefore,  this  test  shows  that  there  is 
spherical  spreading  at  separations  greater  then  4.0  cm. 

Since  the  phase  velocity  is  greater  in  the  sediment-water 
mixture  then  in  water,  the  Rayleigh  distance  is  even  smaller  in  the 
mixture.  These  considerations  establish  that  spherical  spreading  in  the 
sediment  is  a  good  assumption  at  separations  greater  then  4.0  cm. 

3.  Directivity  Since  the  transducers  ere  directional  and  the 
receiving  transducer  is  moved  to  a  new  separation  after  each  set  of 
readings,  greet  cere  has  to  be  taken  so  that  the  transducers  remain  on 
each  others  axes,  where  the  signal  amplitude  is  maximum.  The  axis  is 
the  imaginary  line  running  through  the  center  of  the  circular  element  end 
perpendicular  to  its  face.  The  alignment  problem  is  lessened  if  the 
radiation  pattern  of  the  transducer  is  sufficiently  wide  that  slight 
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errors  in  realigning  the  receiver  exactly  on  the  axis  of  the  transmitter 
will  cause  only  slight  error  in  measuring  the  on-axis  signal  amplitude. 

By  treating  the  transducer  as  a  circular  piston  in  an  infinite, 
rigid  baffle,  the -half-power  angle  can  be  calculated.  If  p(0,r,t)  is  the  on- 
axis  pressure,  then  the  half-power  angle  is  the  angle  at  which  the 
pressure  is  given  by  (Blackstock) 

pfe^.r.t):  p(0,r,t)//2 

For  a  circular  piston,  the  half-power  angle  is  (Blackstock) 

sin  6^  =1.61 6/ka 

In  fresh  water,  at  8  frequency  of  105.0  kHz,  the  calculated  half-power 
angle  for  these  transducers  is  22.6°. 

Figures  4  and  5  show  the  measured  directivity  of  the  two 
transducers  in  fresh  water.  Again,  it  is  seen  that  the  model  of  the 
transducer  as  a  circular  piston  in  a  rigid  baffle  works  well.  The 
directivity  patterns  for  both  transducers  show  that  the  half-power 
angle,  the  angle  at  which  the  pressure  is  3.0  dB  below  the  on-axis 
pressure,  is  about  20°. 

Since  transducer  #1  is  used  as  the  receiver  in  this  experiment, 
and  is  moved  after  each  reading,  it  is  desired  that  the  response  of  this 
transducer  especially  be  constant  for  small  angles  off  axis.  This  is  so 
that  failure  to  exactly  realign  the  transducers  will  still  give  a  reliable 
measurement  of  the  on-exis  value  of  the  signal  amplitude.  Fig.  4  shows 
that  the  sound  pressure  level,  Lp,  for  transducer  #1  varies  less  than 
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0.5  dB  for  7.0*  on  either  side  of  the  axis.  A  0.5  dB  difference  In  Lp 

implies  that  the  actual  pressure  varies  less  than  6.08  from  the  on-axis 
value.  Actually,  the  Lp  for  transducer  *\  is  essentially  constant  for 

angles  up  to  3.0°  off  axis. 

For  media  such  as  saturated  sediments,  in  which  the  sound 
speed  is  higher  than  that  of  water,  the  half-power  angle  would  be  even 
greater.  A  typical  compressions!  wave  speed  for  a  sediment-water 
mixture  is  1800  m/s  (Bell,  1979).  At  the  seme  frequency  that  is  being 
considered,  105.0  kHz,  this  leads  to  a  half-power  angle  of  28*.  Thus  one 
is  led  to  believe  that  for  small  angles  off  axis,  the  signal  amplitude 
would  vary  even  less  in  the  mixture  than  it  would  in  water.  In  the 
measurements  taken  in  this  investigation,  it  was  estimated  that  the 
receiver  was  positioned  in  the  sediment  with  no  more  than  4°  variation 
in  the  angle  off  axis.  Thus,  there  is  confidence  that  the  on-axis  signal 
amplitudes  were  measured  accurately. 

D.  Sediment  Preparation 

The  experimental  sediment  is  made  up  of  MS-MH  size  glass 
beads  manufactured  by  the  Ferro  Corporation  of  Jackson,  Mississippi. 
The  grain  size  distribution  is  shown  in  Fig.  6  (Bell,  1979). 

The  sediment-water  mixture  was  prepared  in  a  five  gallon, 
plastic  bucket  in  the  following  manner.  The  bucket  is  approximately 
35  cm  high  and  29  cm  in  diameter.  The  liquid  in  the  sediment-water 
mixture  is  circulated  by  a  Masterflex  tubing  pump.  Before  the  mixture 
was  put  in  the  bucket,  the  outflow  tube  was  attached  to  the  bottom 


SIZE -mm  1.0 


0.25 


0.0625 


FIGURE  6 

GRAIN  SIZE  DISTRIBUTION  FOR  MS-MH  SIZE  BEADS 
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of  the  bucket  with  duct  tope.  The  bucket  wos  then  partially  filled  with 
de-ionized  water  and  the  input  hose  was  put  in  the  water.  The  pump  was 
turned  on  so  that  all  of  the  air  was  driven  from  the  tube.  To  retard  the 
growth  of  micro-organisms,  20  ml  of  BIO-CLEAN  III  algeecide  was  added 
to  the  water. 

The  glass  shot  was  then  poured  slowly  into  the  water  until  the 
sediment-water  mixture  was  24  cm  deep.  Pouring  wos  halted  at 
intervals  and  the  mixture  kneaded  to  remove  air  trapped  in  the  beads.  A 
layer  of  water,  approximately  6  cm  thick,  wos  left  on  top  of  the 
sediment  mixture  and  the  input  hose  to  the  pump  was  suspended  in  the 
layer,  os  shown  in  Fig.  7. 

The  bucket  wos  then  placed  inside  a  vacuum  tank  on  top  of  a 
vibrator  to  remove  any  air  remaining  in  the  mixture.  The  tank  wos 
evacuated  and  the  mixture  was  vibrated  in  excess  of  6  hours  The 
mixture  was  left  in  the  vacuum  overnight,  but  the  vibrator  was  turned 
off.  After  the  sediment  wos  removed  from  the  vacuum,  the  transducers 
were  inserted  into  the  mixture  so  that  the  center  of  the  elements  were 
12  cm  below  the  sediment-water  interface. 

To  alter  the  fluid  properties,  glycerine  was  poured  into  the 
liquid  on  top  of  the  mixture  and  the  pump  turned  on  until  the  liquid  in  the 
bucket  wos  well  mixed.  Approximately  10  ml  of  algeecide  wos  added 
while  the  liquid  wos  being  circulated  to  further  retard  the  growth  of 
micro-organisms.  The  specific  gravity  of  the  liquid  wos  measured 
several  times  while  the  liquid  wos  being  circulated.  The  liquid  wos 
considered  well  mixed  after  the  values  of  specific  gravity  had  stabilized. 


IV.  EXPERIMENTAL  PROCEDURE 


A.  Phase  Velocitu 

The  oscilloscope  trace  is  triggered  at  the  beginning  of  the 
transmitted  four  cycle  pulse.  The  delay  time  to  the  fourth  positive  peak 
on  the  received  signal  is  measured  with  a  digital  readout  on  the 
oscilloscope.  This  is  first  done  at  a  transducer  separation  of 
approximately  20  cm,  and  then  at  every  2.0  cm  increment  down  to  a 
separation  of  6  cm.  If  these  delay  times  are  then  plotted  on  a  time 
versus  distance  graph,  with  time  as  the  horizontal  axis,  tho  points 
should  fall  on  a  straight  line  parallel  to  the  characteristic  of  the 
wavelet.  The  slope  of  this  line  is  the  phase  velocity. 

Therefore,  to  find  the  phase  velocity  a  least  squares  routine  is 
used  to  fit  the  best  straight  line  through  these  points.  Since  the  line  is 
parallel  to  a  characteristic,  the  equation  of  the  line  is  of  the  form 

r  =  c0t  +  B 

where  c0  is  the  phase  velocity,  B  is  the  intercept  of  the  r-axis,  and  r  and 
t  represent  the  separation  and  time  delay,  respectively.  If  the  time  were 
measured  from  the  fourth  positive  peak  on  the  transmitted  pulse  to  the 
same  peak  on  the  received  pulse,  B  would  ideally  be  zero  (for  a 
nondispersive  medium)  and  the  least  squares  line  would  be  the  same  as 
the  characteristic.  But  since  time  is  measured  from  the  beginning  of  the 
triggered  pulse  and  there  are  additional  time  delays  associated  with  the 
circuitry  and  transducers,  B  turns  out  to  be  negative. 
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This  method  was  tested  using  two  different  size  heads  with 
fresh  water  as  the  pore  fluid:  the  MS-MH  size  which  was  used  in  the 
experiment,  and  the  MS-L  size  which  are  smaller  size  beads.  The  results 
of  these  preliminary  tests  are  shown  in  Figs.  8  end  9.  The  linear 
correlation  coefficient  displayed  in  the  figures  is  a  number  between  zero 
and  one  that  gives  an  indication  of  how  well  the  data  points  fit  a 
straight  line.  A  coefficient  of  one  represents  perfect  linear  correlation 
and  a  coefficient  of  zero  represents  no  correlation.  The  coefficients 
shown  in  the  figures  indicate  a  very  good  correlation  between  the  data 
and  the  straight  lines.  The  good  correlation  suggests  that  the 
disturbances  in  the  sediment,  caused  by  moving  the  transducer  to  each 
separation,  have  a  negligible  effect  on  the  measurement 

In  order  to  check  the  validity  of  our  procedure,  we  measured  the 
phase  velocity  in  water  at  room  temperature,  approximately  25°C,  six 
times.  Each  of  these  results  were  between  1492  m/s  and  1508  m/s.  The 
average  of  these  measurements  was  1501  m/s.  By  linear  interpolation 
of  tabulated  values  of  sound  speed  as  a  function  of  temperature  (AIP 
Handbook,  1963),  the  speed  of  sound  at  25°C  in  fresh  water  is  1496  m/s. 
The  overage  then  has  an  error  of  only  0.3355,  while  the  extreme  values 
have  errors  of  -0.3%  and  0.8%. 

In  another  check,  we  compared  our  results  to  results  obtained  by 
Bell  (Bell,  1979),  Shirley  (Bedford,  Appendix  A,  1982),  and  Hovem  and 
Ingram  (Hovem,  1979).  We  measured  the  phase  velocity  three  separate 
times  in  MS-MH  beads  with  fresh  water  as  the  pore  fluid.  The  three 
results  are  shown  in  Table  I,  along  with  the  other  results  mentioned. 
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BEAD  TYPE:  MS-MH 

FREQUENCY:  105  kHz  VISCOSITY:  0.893  cP 


TIME -ms  SEPARATION -cm 


65.88 

77.63 

88.98 

100.12 

112.07 

121.38 

133.23 

142.31 


5.65 

7.65 

9.65 

11.65 

13.65 

15.65 

17.65 

19.65 


VELOCITY:  1820  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  24  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9995 


FIGURE  8 

PHASE  VELOCITY  IN  MS-MH  BEADS 
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BEAD  TYPE:  MS-L 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


TIME  -  ms 

SEPARATION 

63.62 

5.65 

73.39 

7.65 

83.49 

9.65 

95.51 

11.65 

107.15 

13.65 

118.60 

15.65 

131.73 

17.65 

138.93 

19.65 

VELOCITY:  1791  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  34  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9989 
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FIGURE  9 

PHASE  VELOCITY  IN  MS-L  BEADS 
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Comparison  with  these  other  results  indicates  that  our  results  are  valid. 
But  the  comparison  also  suggests  that  the  velocity  is  strongly  dependent 
on  the  depth  of  the  transducers  in  the  sediment. 

The  comparison  in  the  MS-L  size  beads  is  not  nearly  so  good. 
Our  first  trial  measurement  in  sediments  was  mode  in  the  MS-L  size 
beads.  Since  the  standard  deviation  accompanying  this  measurement  is 
so  large,  34  m/s,  and  since  in  all  subsequent  measurements  the  standard 
deviation  was  substantially  reduced,  the  discrepancy  is  attributed  to 
experimental  error  and  carelessness  in  precisely  determining  the 
separation.  As  the  experiment  progressed,  the  technique  improved  and 
the  errcr  decreased. 


Table  I:  Comparison  of  Velocity  Results  in  MS-MH  and  MS-L  Size  Beads. 


Frequency  -  kHz 

Depth  -  cm 

Velocity  -  m/s 
MS-MH  MS-L 

Bell 

120. 

9 

1809  1854 

Shirley 

1 14. 

10 

1808 

Hovem  &  Ingram 

100. 

20 

1922 

Present  Study 

105. 

12 

1820  1791 

105, 

12 

1840 

105. 

12 

1866 
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Obviously,  the  sound  speed  measurements  in  sediments  are  not 
nearly  as  repeatable  as  the  measurements  made  in  water.  One  reason  for 
this  is  the  uncertainty  in  the  separation  measurement.  In  the  sediment, 
this  uncertainty  is  estimated  to  be  0.1  cm.  At  a  separation  of  20  cm  this 
represents  a  relative  uncertainty  of  0.5  8,  while  at  6  cm  the  uncertainty 
is  1.78.  The  uncertainty  in  the  time  delay  measurement  is  negligible 
compared  to  this.  Relative  uncertainties  of  0.58  and  1.78  in  the 
separation  represent  absolute  uncertainties  of  10  m/s  and  30  m/s  in  the 
velocity,  respectively. 

The  standard  deviations  shown  in  Figs.  8  end  9  were  derived 
from  the  least  squares  method.  In  this  calculation  of  the  uncertainty  it 
was  assumed  that  the  uncertainties  in  the  time  delay  measurements, 
which  were  measured  along  the  horizontal  axis,  were  negligible.  It  was 
also  assumed  that  the  absolute  uncertainties  of  each  of  the  separation 
measurements,  which  were  represented  on  the  vertical  axis,  were  equal. 
Since  both  of  these  assumptions  are  appropriate  and  the  standard 
deviations  calculated  from  the  least  squares  method  are  consistent  with 
uncertainties  in  the  separation  measurements,  these  standard  deviations 
are  a  realistic  measure  of  the  uncertainty  in  the  velocity. 

All  of  the  velocity  measurements  made  in  the  sediments  ere 
listed  in  Table  VII  of  Chapter  IV  end  in  Appendix  0  together  with  the  raw 
data  It  can  be  seen  that  the  standard  deviations  of  the  velocities  were 
less  than  15  m/s  after  the  first  few  trial  measurements.  This  indicates 
that  as  the  experiment  progressed,  the  technique  was  improved  end  the 
uncertainty  decreased. 
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B.  Attenuation 

1.  Measurement  Technique.  In  the  sediment  the  losses  are  due  to 
both  geometrical  spreading  and  absorption  of  the  medium.  Therefore,  to 
calculate  the  attenuation  by  the  medium  the  geometrical  spreading  must 
be  taken  into  account.  As  discussed  in  Section  III.C.2.,  the  transducers 
were  designed  so  that  the  Rayleigh  distance  would  be  relatively  small 
compared  to  the  sediment  sample.  At  distances  greater  than  the 
Rayleigh  distance,  the  signal  amplitude  decays  as  1/r  due  to  spherical 
spreading,  where  r  is  the  separation  between  the  transducers.  This 
distance  was  calculated  for  both  water  and  the  sediment.  Tne 
calculation  for  water  was  then  verified  by  acoustic  measurements. 
These  measurements  support  the  calculated  assumption  that  the 
geometrical  spreading  in  the  sediment  is  indeed  spherical  at  the 
distances  at  which  the  attenuation  measurements  are  token.  Using 
spherical  spreading,  the  losses  due  to  absorption  in  the  sediment  can  be 
determined.  A  detailed  discussion  of  this  procedure  follows. 

The  measurements  establishing  the  1/r  decay  of  the  signal 
amplitude  due  to  spherical  spreading  have  been  discussed  in 
Section  III.C.2.  We  determine  the  absorption  due  to  the  medium  by 
measuring  the  amplitude  of  the  signal  end  assuming  that  the  amplitude 
varies  as 

E  =  (K/r)  e"ar  (IV- 1) 

in  an  absorbing  medium.  Here,  K  is  a  calibration  constant  with  units  of 
voltsxlength  find  <x  is  the  attenuation  in  nepers  per  unit  length.  The  1/r 
factor  in  Eq.  (IV- 1)  takes  into  account  the  spherical  spreading  and  the 
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e'0*1"  term  takes  into  account  the  absorption  of  the  medium.  It  is  also 
assumed  that  the  coupling  between  the  sediment  and  the  transducer  is 
the  same  at  each  separation.  Thus,  the  decay  in  amplitude  of  the  signal 
is  attributed  only  to  the  increase  in  separation. 

The  attenuation  is  determined  by  comparing  the  amplitudes 
measured  at  various  separations.  To  do  this,  one  of  the  amplitude 
measurements  and  the  separation  at  which  it  is  measured  are  designated 
as  reference  values,  Er  and  rr  ,  respectively.  In  this  experiment,  rr  is 

the  smallest  separation  at  which  measurements  are  taken  and  Er  =  E  8t 

r  =  rr  .  Therefore  the  following  relation  holds: 

Er  =  (  K/rr )  e"arr 

From  this  and  Eq.  (IV.  1)  it  follows  that . 

K  =  E  r  ear  =  Er  rr  earr 
eoc(r-rr)  -  Errr/(Er) 

oc(r-rr)  =  In  [Er  rr/(E  r)  1  .  (IV-2) 

Therefore,  if  ln(Errr/(E  r))  is  plotted  versus  (r-rr),  the  points  should  fall 

on  a  straight  line  and  the  slope  of  this  line  is  the  attenuation  To  find  the 
slope  of  this  line,  and  thereby  the  attenuation,  the  method  of  least 
squares  can  again  be  employed. 

i 
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This  technique  was  tested  using  two  different  size  beads  with 
fresh  water  as  the  pore  fluid:  the  MS-MH  size  which  was  used  in  the 
experiment,  and  the  MS-L  size  which  are  smaller  size  beads.  The  results 
of  these  preliminary  tests  are  shown  in  Figs.  10  and  1 1.  As  in  the  case 
of  the  phase  velocity,  a  linear  correlation  coefficient  near  one  indicates 
a  very  close  correlation  between  the  data  and  the  straight  line. 

The  standard  deviation  can  be  calculated  for  this  method, 
assuming  that  the  uncertainties  in  each  of  the  dote  points  ore  equal  and 
that  there  are  no  errors  in  the  values  along  the  horizontal  axis.  Because 
of  these  assumptions,  the  standard  deviation  calculated  from  the  least 
squares  method  is  too  low.  But  the  attenuation  can  be  calculated  by 
other  means,  and  from  these  calculations  a  realistic  estimate  of  the 
error  con  be  made.  These  methods  and  the  uncertainties  associated  with 
them  will  be  discussed  in  the  remainder  of  this  section. 

A  physical  interpretation  of  the  sources  of  error  con  be  acquired 
by  calculating  the  attenuation,  «,  directly  using  Eq.  (IV-2)  with  the  data 
at  any  two  separations.  If  the  data  for  MS-MH  beads  corresponding  to 
separations  of  5.65  cm  and  19.65  cm,  shown  in  Fig.  11,  are  used,  the 
attenuation  calculated  from  Eq.  (IV-2)  is  2.28  Np/m.  We  can  calculate  a 
Pythagorean  estimate  of  the  uncertainty  in  the  attenuation  from  this 
calculation  by  using  the  relation  (Bevington,  Eq  4-9,  1969) 

tfj2  (r-rr)2  =  <?r2(a2+2K/r+ 1 /r2)ref  +  tfr2(a2+2oc/r+ 1  /r2) 

*  (°e/E>2  ref 4  (°e/E)2  (IV-3) 

where  <?E  and  crf  denote  the  uncertainties,  or  standard  deviations,  in  the 


BEAD  TYPE:  MS-L 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893 cP 


AMPLITUDE  -  V 

0.2800 

0.1970 

0.1430 

0.1060 

0.0820 

0.0652 

0.0516 

0.0430 


SEPARATION  -  cm 

5.65 

7.65 

9.65 

11.65 

13.65 

15.65 

17.65 

19.65 


ATTENUATION:  4.71  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.14  Np/m 
LINEAR  CORRELATION  COEFFICIENT:  0.9974 


E- AMPLITUDE  r- SEPARATION 


FIGURE  10 

ATTENUATION  IN  MS-L  BEADS 
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BEAD  TYPE:  MS-MH 


FREQUENCY:  105.0  kHz 

VISCOSITY:  0.893 cP 

AMPLITUDE  -  V 

SEPARATION  -  cm 

0.3720 

5.65 

0.2690 

7.65 

0.2004 

9.65 

0.1570 

11.65 

0.1300 

13.65 

0.1074 

15.65 

0.0692 

17.65 

0.0777 

19.65 

ATTENUATION: 

2.40  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.09  Np/m 
LINEAR  CORRELATION  COEFFICIENT:  0.9962 


FIGURE  11 

ATTENUATION  IN  MS-MH  BEADS 
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AS-85-63 

RDC-GA 

1-23-85 


Amplitude  ond  separation,  respectively,  ond  the  quantities  mocked  "ref" 
signify  thot  the  seporotions  ond  Amplitudes  correspond  to  the  seporotion 
of  5.65  cm.  To  perform  this  calculation,  the  uncertainties  in  the 

separation  and  the  amplitude,  or  and  o^,  must  first  be  determined.  The 

uncertainty  in  the  separation  is  estimated  to  be  0.1  cm. 

The  uncertainty  in  the  amplitude  has  two  sources.  The  first  one 
is  instrumental,  because  the  oscilloscope  can  be  read  accurately  to  1/4 
of  a  gradation.  The  scale  is  changed  at  different  separations,  since  the 
amplitude  grows  with  decreasing  r,  so  this  type  of  uncertainty  increases 
with  decreasing  separation. 

The  second  source  of  uncertainty  in  E  is  due  to  the  uncertainty 
in  r.  Since  the  amplitude  is  a  function  of  r,  as  was  assumed  in  Eq.  (IV- 1), 
any  uncertainty  in  r  produces  an  uncertainty  in  E.  This  type  of  error  can 
be  estimated  from  the  relation  (Bevingtori,  Eq.  4-9,  1969) 

~  °r  dE/dr  =  ar  E  (  ar  +  1j7  J 

By  calculating  this  uncertainty,  and  comparing  the  results  with  the 
instrumental  uncertainty  in  E,  it  can  be  seen  that  over  60%  of  the 
uncertainty  in  E  is  due  to  the  uncertainty  in  r.  Using  Eq.  (IV-3),  and  the 
values  in  the  rows  in  Table  II  corresponding  to  r=  19.65  cm  and  r=5.65  cm, 
the  uncertainty  in  the  direct  calculation  is  found  to  be  0.522  Np/m.  Note 

from  Eq.  (3)  that  the  uncertainty,  Oj,  is  inversely  proportional  to  the 

difference  between  the  two  separations  being  compared.  Thus,  one  way 
to  improve  the  measurements  would  be  to  increase  the  distances  over 
which  the  measurements  are  taken. 


Using  each  pair  of  data  shown  in  Fig.  1  1,  the  attenuation  can  be 
calculated  seven  times  in  this  way,  and  the  values  designated  as 
These  seven  values  of  attenuation  can  then  be  weighted  with  their 
corresponding  uncertainties,  and  averaged  as  shown  in  Table  II 

(Bevington,  p.  73,  1969).  The  resulting  average  is  a  better  estimate  of 
the  attenuation  than  any  of  the  individual  attenuations,  and  the  standard 
deviation  calculated  from  this  method  is  a  better  estimate  of  the 
uncertainty.  Although  this  value  of  a  differs  by  5%  from  the  least 
squares  value,  this  standard  deviation  is  a  more  realistic  estimate  of  the 
uncertainty. 


Table  II:  Analysis  of  Errors  in  Attenuation. 


r  -  cm 

E  -  volts 

oE/e 

<Jr/(0C2+2(X/r+  Vf2) 

cr1  -  Np/m 

-Np/m 

5.65 

0.372 

0.0380 

0.000404 

7.65 

0.269 

0.0407 

0.000239 

3.06 

1.057 

9.65 

0.2004 

0.0460 

0.000163 

1.606 

2.082 

11.65 

0.157 

0.0429 

0.000121 

1.029 

2.317 

13.65 

0.130 

0.0476 

0.000095 

0.813 

2.116 

15.65 

0.1074 

0.0533 

0.000077 

0.690 

2.235 

17.65 

0.0892 

0.0592 

0.000065 

0.613 

2.408 

19.65 

0.0777 

0.0587 

0.000056 

0.522 

2.283 

1 

<?K2  =  -  so  that  =  0.296  Np/m 

2  ('/eg2 


oc  =  <?K2  ]>  (ai/<7j2)  =  2.266  Np/m 


In  calculating  the  attenuation  using  this  method,  the  error  in  the 
reference  separation  and  amplitude  is  repeated  in  each  of  the  seven 
calculations.  Thus,  if  reference  measurements  are  particularly  bad,  the 
error  in  them  will  be  compounded  in  the  final  result.  To  avoid  this  type 
of  error,  the  attenuation  can  be  calculated  with  any  two  pairs  of  data,  E 
and  r.  The  values  corresponding  to  r  =  5.65  cm  do  not  necessarily  need  to 
be  the  reference  separation  and  amplitude.  Since  there  are  eight  pairs  of 
data,  there  are  28  possible  ways  to  calculate  <x  using  Eq.  (IV-2).  The 
unweighted  mean  of  these  28  attenuations  and  the  standard  deviation  of 
this  mean  can  be  easily  calculated.  For  the  MS-MH  beads  of  Fig.  1 1,  this 
mean  is  «  =  2.38  Np/m  and  the  standard  deviation  is  crK  =  0,47  Np/m. 

This  value  of  the  attenuation  differs  only  by  0.6*  from  the  least  squares 
value.  But  this  standard  deviation  is  larger  than  the  one  calculated  from 
the  weighted  averages.  This  may  signify  that  the  estimate  of  the 

uncertainty  in  r  used  in  that  calculation,  <?r  =  0.1  cm,  was  too  low.  In 

later  calculations,  however,  there  is  better  agreement  between  these 
two  calculations  of  the  standard  deviation.  This  suggests  that  the 
separation  was  determined  more  precisely  as  the  experiment  progressed. 
Another  reason  for  the  discrepancy  is  that  all  of  the  attenuations  are 
weighted  equally  in  calculating  the  mean.  But  as  can  be  seen  from 
Eq.  (IV-3),  the  uncertainty  decreases  as  the  difference  between  the  two 
separations  increases.  Thus,  the  standard  deviation  calculated  by  this 
method  represents  an  upper  bound  of  the  actual  uncertainty  in  the 
attenuation. 
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The  attenuation  was  calculated  for  each  value  of  viscosity  using 
the  method  of  least  squares  and  the  mean  of  all  combinations.  In  each 
case,  the  two  values  of  attenuation  differed  by  less  than  18.  But  we  feel 
that  the  standard  deviation  associated  with  the  mean  of  ell  combinations 
was  more  realistic. 

It  can  be  noted  that  the  uncertainty  in  the  separation  is  by  far 
the  biggest  source  of  error.  Therefore,  to  improve  these  measurements, 
a  method  for  determining  the  separation  of  the  transducers  more 
accurately  would  have  to  be  devised.  Toward  this  end,  the  separations 
were  corrected  using  the  velocity  results.  The  attenuations  were  then 
calculated  by  the  mean  of  all  combinations  with  these  corrected 
separations.  These  attenuations  were  within  18  of  the  attenuations 
calculated  by  the  other  methods,  except  for  the  case  of  the  NS-MH  beads 
of  Pig.  1 1  where  the  difference  was  2.58.  In  all  but  a  few  cases,  the 
standard  deviation  was  substantially  reduced.  Thus,  the  attenuation  was 
calculated  by  the  mean  of  all  combinations.  In  all  but  those  few  cases 
where  the  standard  deviation  was  not  reduced,  the  separations  were 
corrected  using  the  velocity  results.  Reducing  the  data  in  this  way,  we 
are  confident  that  the  attenuations  were  determined  within  3  dB/m, 
except  for  a  few  cases.  The  attenuations  and  their  uncertainties, 
corresponding  to  each  fluid  viscosity,  ere  listed  in  Table  VII  in 
Chapter  IV.  They  are  also  listed  together  with  the  raw  data  in 
Appendix  B. 
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2.  Comparison  of  Results.  To  determine  the  validity  of  the 
attenuation  results,  the  measurements  in  the  water  saturated  beads, 
with  no  glycerine  added,  were  compared  to  the  results  of  David  Bell 
(Bell,  1979).  His  attenuation  measurements  were  given  in  decibels  per 
wavelength.  Converting  his  results  to  nepers  per  meter  at  105.0  kHz 
resulted  in  the  values  of  5.45  Np/m  and  3.13  Np/m  for  MS-L  and  MS-MH 
size  beads,  respectively. 

The  discrepancy  in  attenuation  between  Bell's  values  and  those 
presented  in  Figs.  10  and  11  was  0.74  Np/m  in  the  first  case  and 
0.73  Np/m  in  the  second.  Since  this  difference  was  too  large  to  be 
attributed  to  experimental  error,  an  effort  was  made  to  explain  it. 

In  order  to  account  for  the  geometrical  spreading  of  the 
transducers,  Bell  measured  the  amplitude  decay  in  water  at  120  kHz  and 
the  amplitude  decay  in  the  sediment  at  the  some  frequency.  He  assumed 
that  the  only  losses  in  the  water  were  due  to  geometrical  spreading  and 
that  the  losses  in  the  sediment  were  due  to  geometrical  spreading  and 
absorption  in  the  medium.  He  then  assumed  that  the  losses  due  to 
geometrical  spreading  in  the  sediment  were  the  same  as  in  water,  at  the 
same  frequency.  Therefore,  to  calculate  the  attenuation  he  subtracted 
the  spreading  losses,  which  he  determined  from  his  measurements  in 
water,  from  the  total  losses  (Bell,  1979). 

In  doing  this,  he  assumed  that  the  geometrical  spreading  in  the 
water  was  the  same  as  in  the  sediment,  for  measurements  taken  at  the 
some  frequency.  But  the  Rayleigh  distance,  which  is  a  measure  of  the 
distance  at  which  spherical  spreading  starts,  is  inversely  proportional  to 
the  sound  speed  Since  the  sound  speed  is  greater  in  the  sediment  than  in 


the  water,  the  Rayleigh  distance  is  smaller  in  the  sediment  than  in  the 
water.  Also,  the  distances  at  which  he  took  his  measurements  were  less 
than  the  Rayleigh  distance  for  his  transducers,  so  that  in  water  his 
amplitude  did  not  vary  as  1/r  Therefore,  the  losses  due  to  geometrical 
spreading  were  less  in  the  water  than  in  the  sediment.  The  extra 
geometrical  spreading  that  was  occurring  in  the  sediment  was  attributed 
to  absorption  in  the  sediment,  resulting  in  a  value  of  attenuation  that 
was  too  large. 

We  also  compared  our  measurements  in  MS-MH  beads  with  those 
made  by  Shirley  (Bedford,  Appendix  A,  1982)  and  Hovem  and  Ingram 
(Hovem,  1979).  These  results  are  summarized  in  Table  III.  Shirley's 
result,  like  Bell's,  is  much  higher  than  ours  We  believe  that  he  measured 
the  attenuation  in  the  same  way  that  Bell  did.  Thus,  we  feel  that  his 
results  differ  from  ours  for  the  same  reason  that  Bell's  did. 

Hovem  and  Ingram's  results  at  100  kHz  also  differ  significantly 
from  ours.  If  their  results  at  just  this  one  frequency  are  considered,  one 
might  be  led  to  the  conclusion  that  the  measured  attenuation  is  very 
much  dependent  on  the  depth  of  the  transducers  in  the  sediment.  But  if 
their  results  at  other  frequencies  in  that  range  ere  considered,  it  can  be 
seen  that  their  data  exhibit  a  significant  amount  of  frequency 
dependence.  It  is  difficult  to  say  whether  this  exhibited  frequency 
dependence  is  due  to  actual  physical  processes  or  to  scatter  in  the  data. 
Other  frequency  dependent  measurements  need  to  be  taken  to  investigate 


Table  III:  Comparison  of  Attenuation  Results  in  liS-MH  Size  Beads. 


Frequency  -  kHz 

Depth  -  cm 

Attenuation 

Present  Study 

105.0 

12. 

20.85 

20.74 

21.61 

Bell 

120.0 

9. 

27.10 

Shirley 

1 14.0 

10. 

30.91 

Hovem  & 

60.0 

20. 

24.86 

Ingram 

90.0 

20. 

17.63 

100.0 

20. 

13.20 

150.0 

20. 

19.65 

There  were  several  assumptions  made  in  our  measurements 
which  could  affect  the  validity  of  the  attenuation  results.  One 
assumption  was  that  the  mechanical  coupling  between  the  transducer  and 
the  sediment  was  the  same  at  each  separation.  To  justify  this 
assumption,  each  time  the  transducer  was  moved  the  sediment-water 
mixture  was  vibrated  gently  for  approximately  thirty  minutes.  The 
mixture  was  then  left  undisturbed  for  another  half  hour  before  another 
measurement  was  taken.  Four  or  five  more  readings  were  then  taken 
over  the  next  three  to  four  hours.  Since  we  found  that  the  amplitude  and 
time  delay  readings  did  not  change  over  this  period  of  time,  except  for 
small  variations  attributed  to  experimental  uncertainty,  it  was  assumed 
that  the  coupling  between  the  transducer  and  sediment  did  not  vary. 


3.  Transient  Considerations.  From  a  theoretical  viewpoint,  it 
would  be  desirable  to  measure  the  time  delay  and  amplitude  using  a 
steady  state  signal.  In  this  way,  we  could  be  confident  that  the  signal 
had  a  single  frequency.  The  transient  characteristics  of  the  transducers 
could  also  be  neglected.  However,  the  container  holding  the  sediment  has 
a  limited  volume.  Thus,  there  will  be  reflections  of  the  signal  from  the 
sides  of  the  container  and  from  the  sediment-water  interface.  To 
prevent  these  reflections  from  interfering  with  the  signal  that  travels 
along  the  direct  path  between  the  two  transducers,  a  four  cycle  pulse 
was  transmitted.  The  signal  was  then  cut  off  for  several  milliseconds 
so  that  all  the  reflections  in  the  sediment  would  attenuate  before  a  new 
pulse  was  transmitted.  A  photograph  of  the  oscilloscope  trace  of  the 
transmitted  and  received  signals  is  shown  in  Fig.  12.  Note  that  the 
received  signal  contains  more  than  four  cycles.  Since  the  transducers 
are  being  operated  at  one  of  their  resonant  frequencies,  they  continue  to 
vibrate,  or  ring,  for  several  cycles  after  the  input  signal  is  turned  off. 

In  this  section  and  the  next,  an  attempt  is  made  to  justify 
treating  the  pulses  as  steady  state,  harmonic  signals  in  analyzing  the 
data  and  comparing  it  to  the  theory.  This  section  discusses  the  transient 
effects  of  the  transducers  on  the  signal  amplitude.  The  signal  amplitude 
can  also  be  effected  by  dispersion.  These  effects  are  discussed  in  the 
next  section. 

Since  we  compared  amplitude  measurements  at  different 
separations  to  calculate  the  attenuation,  it  is  assumed  that  the 
measured  voltage  amplitude  is  linearly  proportional  to  the  acoustic 
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HORIZONTAL  SCALE:  20*is/DIV 

VERTICAL  SCALE: 

(a)  RECEIVED  SIGNAL -0.05  V/DIV 

(b)  TRANSMITTED  SIGNAL -5  V/DIV 

SEPARATION:  8.04cm 


FIGURE  12 

TRANSMITTING  AND  RECEIVING  PULSES  IN  SEDIMENTS 
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pressure  amplitude  that  produced  it  in  the  sediment.  In  order  to  justify 
this  assumption,  consider  the  block  diagram  in  Fig.  13.  The  amplitude  of 
the  vollage  signal  put  in  to  the  transmitter  is  denoted  by  Ej.  This 

produces  an  acceleration  of  the  transducer  face,  denoted  by  u.  The  signal 
produced  by  this  motion  of  the  transducer  face  travels  a  distance  r  to  the 
receiver,  where  it  produces  a  pressure  P0.  The  receiver  converts  this 
pressure  signal  to  an  output  voltage  signal  E0.  The  phase  constant  <f> 
signifies  that  the  measurements  are  taken  on  the  fourth  cycle.  Thus 
E0(f),  shown  in  Fig.  12,  represents  the  peak-to-peak,  voltage  amplitude 
of  the  fourth  cycle  of  the  received  signal.  The  operators  Lj,  L2,  arid  L3 

denote  linear  functions  of  their  arguments. 

Our  use  of  the  phase  constant  <f>  deserves  some  explanation.  If  T 
represents  the  period  of  one  cycle,  then  the  positive  and  negative  peaks 
of  the  fourth  cycle  will  correspond  to  the  phases 

+♦  =  U  -  r/Co  *  3.25  T 

and  «j>_  =  t_  -  r/c.Q  =  3.75  T  , 

respectively.  Note  that  the  separation  r  is  the  same  in  each  case,  so  that 
the  difference  in  the  travel  times  t+  and  t_  is  0.5T.  The  pressure 
amplitude  of  the  fourth  cycle  is  then 

P„(t)  =  p0(<h>  -  po<i> 

Since  u  represents  the  acceleration  of  the  transducer  face,  and  the  phase 
represented  by  spans  half  of  a  period,  the  quantity  u(<(>)  does  not  really 
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FIGURE  13 

LINEARITY  OF  THE  TRANSDUCERS 


make  sense.  But  in  Fig.  13  it  represents  the  quantity 

u(<f>)  =  u(<(>+)  -  u(<j> J 

The  terms  £0(<f>)  and  have  similar  meanings. 

Consider  the  two  graphs  of  E0(f)  versus  E|(<j>)  shown  in  Fig.  14. 

These  were  taken  by  varying  the  amplitude  of  the  input  voltage,  Ej(<j>), 

and  measuring  the  amplitude  of  the  fourth  cycle  of  the  output  voltage 
signal.  The  graphs  were  made  at  the  distances  r  =  10.0  cm  and 
r  =  20.0  cm  so  that  the  entire  range  of  received  amplitudes  that  were 
measured  in  the  experiment  would  be  covered.  The  graphs  show  that  the 
output  amplitude  of  the  fourth  cycle,  E0(<j>),  is  a  linear  function  of  the 

input  amplitude  Ej.  This  can  be  expressed  as 

E0<+)  =  L3(Et(4>)) 

To  calculate  the  attenuation,  it  is  assumed  that  the  pressure 
signal  decays  according  to 

P0(r,t)  =  (K/r)  e'"Kr  u(t-r/Co) 

where  K  represents  an  arbitrary  constant,  and  c0  is  the  phase  velocity. 
The  axial,  farfield  pressure  due  to  a  transient  signal  in  a  nonabsorbing, 
nondispersive  medium  is 

P0(r,t)  =  (K/r)  u(t-r/Co)  • 

To  account  for  the  attenuation,  the  exponential  term  was  incorporated, 
ad  hoc.  Taking  the  attenuation  into  account  in  this  way  is  justified 
because  the  farfield  pressure  on  the  axis  for  steady  state,  harmonic 
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(a)  LINEAR  AMPLITUDE  RELATIONSHIP  AT  20  cm 


(b)  LINEAR  AMPLITUDE  RELATIONSHIP  AT  10 cm 


FIGURE  14 

INPUT  VOLTAGE  AMPLITUDE  versu*  OUTPUT  VOLTAGE 
AMPLITUDE  OF  THE  FOURTH  CYCLE 
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radiation  from  a  baffled  piston  is  given  by  (Blackstock) 

P0(r,t)  =  (K/r)  ei(u)t“kr) 

The  e“ar  term  is  arrived  at  in  this  case  by  assuming  that  the  wave 
number  k  is  complex,  where  a  * s  its  imaginary  part.  Since  u(<j>)  is  a 
constant,  we  get 

P0($)  =  (K/r)  e"ar 

Therefore,  the  assumption  used  to  determine  the  attenuation, 
E0(4>)  =  K/r  e-«r  ,  (1) 

will  be  valid  if  E0(^)  =  l2(Po($))  •  (2) 

To  show  that  (1)  is  valid,  assume  that  (2)  is  true  and  also  that 

u(f)  =  Lf(Ef<f» 

This  is  the  same  as  assuming  that  the  transducers  produce  a  signal 
amplitude  that  is  linearly  proportional  to  their  signal  input,  for  a  given 
cycle,  it  therefore  follows  that: 

E„($) « l2(P0<t» 

=  L2  ((K/r)  e_c<r  u(f )] 

=  (K/r)  e“ar  L2(u(f )) 

=  (K/r)  e-«r  L2{L1[E1(<)>))} 

=  (K/r)  e_c<r  L3lE1<4>)] 

=  (K'/r)  e"ar  Et(<t>) 

where  K'  is  another  arbitrary  constant.  For  a  fixed  r,  this  shows  that 
E0(<j>)  is  linearly  proportional  to  Ej(<|>),  which  has  been  demonstrated  in 
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the  two  graphs  of  Fig.  14.  This  result  supports  the  assumption  (2). 
Furthermore,  Ej(<j>)  is  a  constant  so  that  (1)  directly  follows. 

4.  Effects  of  Dispersion.  In  determining  the  attenuation,  the 
decrease  in  amplitude  with  increasing  separation  was  attributed  to  two 
causes;  spherical  spreading  and  absorption  of  the  medium.  In  reducing 
the  data,  it  was  assumed  that  the  signal  contained  only  one  frequency 
and  that  any  dispersion  could  be  neglected.  But,  since  the  transmitted 
signal  was  the  product  of  a  sine  wave  and  the  difference  between  two 
step  functions,  the  signal  contained  a  continuum  of  frequencies.  We 
justified  treating  the  signal  as  a  pure  tone  by  channeling  the  signal 
through  a  narrow  bandpass  filter  and  by  making  our  measurements  on  the 
fourth  cycle  of  the  signal  to  minimize  transient  effects. 

Since  the  signal  did  contain  more  than  one  frequency,  dispersion 
would  cause  a  decrease  in  amplitude  with  increasing  separation. 
However,  we  believe  that  the  effect  of  dispersion  on  our  measurements 
is  very  small.  The  theory  and  the  measurements  predict  little  change  in 
the  phase  velocity  of  the  compressional  wave  of  the  first  type  with 
frequency.  Since  the  received  signal  is  filtered  as  a  further  precaution, 
we  would  expect  the  influence  of  dispersion  on  the  amplitude  to  be 
negligible. 

To  be  more  confident  in  this  assumption,  the  attenuation  was 
determined  by  measuring  the  amplitudes  from  three  consecutive  cycles 
on  the  received  signal:  the  third,  fourth,  and  fifth  cycles.  This  was  done 
for  three  different  sets  of  measurements.  The  results  are  shown  in 
Table  IV.  In  each  set,  the  three  attenuations  are  within  a  standard 


deviation.  The  attenuation  of  the  middle  set  exhibits  the  most  spread, 
but  even  so  the  spread  is  still  within  the  estimated  uncertainties. 
Curiously,  the  standard  deviations  corresponding  to  this  set  are  the 
smallest  in  the  table. 

If  dispersion  was  a  significant  factor,  the  envelope  of  the 
received  signal  would  change  shape  as  the  signal  propagated  through  the 
sediment.  This  would  cause  the  attenuations  measured  from  amplitudes 
of  different  cycles  to  be  distinctly  different.  Since  they  ere  not,  we 
conclude  that  any  error  in  the  attenuation  caused  by  dispersion  is 
negligible  within  the  limits  of  error  of  our  measuring  technique. 


Table  IV:  Comparison  of  Attenuations  and  Standard  Deviations 

Measured  from  Different  Cycles  of  the  Same  Received  Signal. 


CYCLE 


3 

4 

5 


3 

4 

5 


3 

4 

5 


VISCOSITY  -cP 


1.33 


1.33 


2.59 


ATTEN  -  dB/m  STD  DEV  -  dB/m 


20.5 


20.4 

19.5 


20.7 

19.0 

18.1 


23.6 

23.2 

23.5 
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2.4 


■  -A.-  t.m 


V.  EXPERIMENTAL  RESULTS 


A.  Phusical  Properties  of  the  Sediment 

In  order  to  calculate  the  wave  number  using  the  theory  outlined 
in  Chapter  II,  the  physical  properties  of  the  saturated  sediment  needed  to 
be  determined.  The  values  of  most  of  these  parameters  were  easily 
obtained,  either  by  direct  measurement  or  from  tables  listing  material 
properties,  which  were  found  in  handbooks  and  journal  articles. 
Assigning  numerical  values  to  some  of  the  other  parameters,  on  the  other 
hand,  was  not  so  straightforward. 

The  physical  properties  of  the  solid  and  fluid  constituents  were 
the  most  easily  determined.  Glass  beads  constitute  the  solid  part  of  the 
mixture.  The  density  and  bulk  modulus  of  glass  were  found  in  several 
sources  (Bedford,  1982;  Kinsler,  1962).  The  average  grain  size,  or  bead 
diometer,  was  measured  by  Bell  from  a  grain  size  analysis.  The  porosity 
of  MS-MH  size  beads  was  also  determined  by  Bell  (Bell,  1979).  These 
values  are  listed  in  Table  V. 

The  properties  of  the  fluid  change  in  the  experiment  since 
glycerine  is  added  to  the  mixture  after  each  set  of  measurements.  The 
density  of  the  fluid  is  easily  determined  by  measuring  the  specific 
gravity  with  a  hydrometer.  Knowing  the  specific  gravity  and  the 
temperature,  the  viscosity  was  read  from  a  table  containing  values  of 
viscosity  for  solutions  of  varying  glycerine-water  concentrations 
(Sheely,  1932).  Shirley  measured  the  sound  speed  in  glycerine-water 


Table  V:  Physical  Properties  and  Parameters 

of  Sediment 

Used  in 

Theoretical  Calculations. 

PARAMETER  SYMBOL 

PRESENT 

STUDY 

BEDFORD 

SHIRLEY 

Porosity 

♦ 

0.365 

0.365 

0.365 

Bead  Diameter  (cm) 

dm 

0.0177 

0.018 

0.0177 

Glass  Bulk 

Modulus  (*10"  dyn/cm2) 

KS 

3.9 

4.07 

3.5 

Glass  Density  (g/cm3) 

Ps 

2.5 

2.5 

2.5 

Frame  Bulk 

Modulus  (dyn/cm2) 

*b 

8.0x1 07 

5.8x1 03 

8.0x1 07 

Frame  Shear 

Modulus  (dyn/cm2) 

^b 

8,0x1 07 

5.8x1 09 

8.0x1 07 

Comp.  Log  Decrement 

5c 

0.1 

0.15 

0.1 

Shear  Log  Decrement 

5s 

0.1 

0.15 

0.1 

Permeability  (xio~?  cm2) 

Bo 

3.8 

2.48 

2.64 

Kozeny-Carman  Constant 

ko 

2.76 

4.23 

3.98 

Depth  in  Sediment  (cm) 

12 

20 

10 
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solutions  with  varying  concentrations  of  glycerine  (Bedford,  Appendix  A, 
1982).  From  this  data,  and  values  of  density  and  viscosity  at  20*C  found 
in  the  CRC  Handbook,  he  was  able  to  determine  the  bulk  modulus  as  a 
function  of  viscosity  for  concentrations  up  to  Z2%  glycerine,  by  weight. 
But  since  his  measurements  were  taken  at  room  temperature,  which  is 
closer  to  25°C,  we  adjusted  his  measurements  so  that  the  bulk  modulus 
values  corresponded  to  viscosity  values  at  25°C.  His  results  were  then 
projected  to  include  glycerine  concentrations  up  to  40£.  The  results  of 
this  procedure  are  shown  in  Table  VI.  This  extrapolation  of  his  date  is 
questionable  at  best,  but  ihere  was  not  enough  time  to  repeat  the 
measurements.  These  adjustments  produce  only  slight  changes  in  the 
theoretical  predictions.  At  some  point  in  the  future  it  would  be 
desirable  to  redo  and  hopefully  repeat  Shirley's  measurements 

The  properties  of  the  solid  matrix  are  more  difficult  to 
determine.  With  all  of  the  other  properties  defined,  these  properties 
were  chosen  so  that  the  theoretical  predictions  and  the  experimental 
data  agreed  for  the  case  of  pure  water.  Once  this  agreement  was 
satisfied,  we  were  able  to  show  that  the  Biot-Stoll  theory  predicted  the 
same  change  in  attenuation  and  phase  velocity  that  was  indicated  by  the 
experimental  data  as  the  viscosity  of  the  pore  fluid  was  increased. 

Four  of  these  properties  were  the  real  and  imaginary  parts  of 
the  bulk  end  shear  moduli  of  the  drained  glass  beads.  Hovem  and  Ingram 
(Hovem,  1979),  Shirley  (Bedford,  Appendix  A,  1982),  and  Bedford 
(Bedford,  1984)  chose  the  real  parts  of  both  of  these  moduli  so  that 
ther°  was  agreement  between  the  theory  and  compressionol  velocity 
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Table  VI:  Physical  Properties  of  Fluid  Used  in  Theoretical  Calculations. 


VISCOSITY 

FLUID  DENSITY 

FLUID  BULK  MODULUS 

T[  -  CP 

pf  -  g/cm3 

Kf  -  *1010  dyn/cm2 

0.893 

0.9971 

2.175 

1.0 

1.0080 

2.284 

1.2 

1.0242 

2.442 

1.4 

1.0373 

2.574 

1.6 

1.0433 

2.683 

1.8 

1.0576 

2.781 

2.0 

1.0654 

2.867 

2.2 

1.0723 

2.943 

2.4 

1.0785 

3.012 

2.6 

1.0840 

3.075 

2.8 

1.0887 

3.128 

date.  Bedford  made  theoretical  comparisons  with  Hovem  and  Ingram’s 
data  and  his  choices  of  the  parameters  ore  very  similar  to  Hovem  and 
Ingram’s. 

Bedford  had  good  agreement  between  theory  and  data  for 
Kb  =  Mb  =  5.6*109dyn/crr|2 

But  Hovem  and  Ingram's  measured  compressional  velocity  at  1 00  kHz, 
1922  m/s,  is  considerably  higher  than  ours  for  the  pure  water  case, 
1830  m/s.  This  is  presumably  because  their  measurements  were  made  8 
cm  deeper  in  the  sediment  than  ours  were.  As  a  result,  their  choice  of  Kb 

and  produced  theoretical  velocities  that  were  too  large  for  our  data. 

But  the  theoretical  attenuations  produced  by  this  choice  compared  well 
with  our  data.  Shirley’s  measured  velocity,  at  1 14  kHz  with  pure  water 
as  the  pore  fluid,  agreed  quite  well  with  our  results.  It  should  be  noted 
that  our  measurements  were  taken  at  similar  depths.  His  choice  for  the 
two  moduli  was 

Kb  =  «b  =  8.0  x  10f  d«'Vcm2 

This  choice  produced  theoretical  velocities  that  agreed  quite  well  with 
our  data.  The  disparity  between  these  two  sets  of  values  suggests  that 
the  properties  of  the  sediment  change  drastically  with  increasing  depth 
in  the  sediment. 

We  used  Shirley’s  values  for  the  bulk  and  shear  moduli,  even 
though  the  theoretical  attenuations  they  produced  were  much  too  high. 
This  discrepancy,  whicn  was  the  result  of  differences  in  experimental 
technique  and  data  reduction,  was  discussed  in  Section  IV  B  2.  In  order 
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to  achieve  a  good  fit  between  the  theory  and  the  experimental  data  for 
the  attenuation  case.,  we  had  to  choose  other  parameters  differently 
from  Shirley. 

The  imaginary  parts  of  the  bulk  and  shear  moduli  of  the  granular 
matrix  are  the  terms  that  take  into  account  the  losses  at  the 
grain- to-grain  contacts.  From  our  experience,  these  parameters  have 
negligible  effect  on  the  compressional  velocity,  and  only  a  slight  effect 
on  the  attenuation.  In  the  literature,  these  two  parameters  have  been 
related  to  their  real  counterparts  by  two  other  parameters,  the 

compressional  and  shear  logarithmic  decrements,  Sc  and  Ss,  respectively. 

The  relationships  defining  these  parameters  are 

Kb'  =  sc  ^b/^  and  Mb'  =  5s  ^b/^ 

where  a  primed  variable  denotes  the  imaginary  part  of  the  modulus  and 
an  unprimed  varible  denotes  the  real  part.  Our  choices  of  these  variables 
are  listed  in  Table  V. 

The  last  parameter  to  be  chosen  was  the  permeability.  Bell 
measured  this  quantity  directly  for  four  sizes  of  glass  beads  and  three 
different  sands  (Bell,  1979).  For  the  MS-MH  size  beads  he  measured  the 
permeability  to  be  B0  =  1.70  *  10"7  cm2 .  The  permeabllty  is  related  to 
an  empirical  parameter,  the  Kozeny-Carmen  constant,  k0,  by  the  relation 
(Bell,  1979) 

V*3 

= - 

36  k0  (1-<)>)2 
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Bell's  value  for  B0  corresponded  to  k0  =  6.17.  Hovem  and  Ingram  chose 
k0=  4.23  so  that  B0  =  2.48  *  10  ~7  cm2  Shirley  used  k0  =  3.98  so  that 
B0  =  2.64  x  10  cm2.  We  found  good  correlation  between  theory  and 
data  by  choosing  B0  =  3.7  x  io  "?  cm2.  This  corresponded  to  k0  =  2.83. 
This  was  within  the  range  of  values  in  the  data  presented  by  Stoll  (Stoll, 
1974)  and  also  within  the  range  of  values  measured  by  Bell  for  several 
sands  and  sizes  of  glass  beads. 

Table  V  summarizes  the  values  of  the  physical  properties  and 
parameters  used  in  the  present  study,  and  also  those  chosen  by  Bedford 
(Bedford,  1984)  and  Shirley  (Bedford,  Appendix  A,  1982).  In  Fig.  15  our 
data  are  compared  with  the  theory  with  Bedford's,  Shirley’s,  and  our 
choices  of  the  parameters. 
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Table  VII:  Velocity  and  Attenuation  in  MS-MH  Beads  with 
Variable  Fluid  Properties. 


VISCOSITY 

VELOCITY 

STD.  DEV. 

ATTENUATION 

STD  DEV. 

T\  -  CP 

c0  -  m/sec 

aC()  -  m/sec 

<x  -  dB/m 

crK  -dB/m 

DATPU  1  __ 

Dn  1  Lril  1 

0.89 

1820.0 

23.6 

20.7 

4.1 

0.89 

1869.7 

9.0 

20.7 

1.1 

1.07 

1861.3 

6.8 

21.3 

2.2 

1  >22 

1878.9 

2.6 

21.0 

3.1 

1.33 

1893.3 

4.2 

20.4 

2.9 

1.33 

1898.6 

5.7 

19.0 

1.9 

—  DATPU  O  - 

Dn  1  Ln  Z. 

0.89 

1839.4 

15.2 

21.7 

2.5 

1.15 

1873.8 

8.7 

23.1 

2.2 

1.48 

1966.7 

9.4 

26.0 

3.6 

1.74 

1973.3 

8.9 

20.9 

2.8 

2.16 

2044.1 

7.3 

21.0 

1.7 

2.59 

2067.5 

9.2 

23.2 

2.1 

2.81 

2074.5 

13.0 

24.3 

1.4 

3.23 

2098.8 

12.1 

24.3 

2.2 

-  DILUTING  BATCH 

2  - 

2.46 

2034.7 

8.2 

26.0 

1.7 

2.09 

1986.6 

3.6 

20.8 

1.6 
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repeatability  of  the  results.  The  results  for  the  three  cases  agree 
reasonably  well  with  each  other  and  also  with  the  theoretical  curves. 
The  figures  demonstrate  that  these  experimental  results  are  repeatable 
within  the  estimated  uncertainties,  especially  in  the  phase  velocity 
case.  The  larger  amount  of  scatter,  and  the  larger  uncertainties,  in  the 
attenuation  results  suggest  that  further  improvements  can  be  made  in 
that  measurement  technique.  One  suggestion  would  be  to  design  a 
bracket  that  would  hold  the  transducer  mount  and  guide  it  straight  into 
the  sediment  while  it  was  being  repositioned  This  would  probably  be 
the  easiest  and  least  expensive  way  to  reduce  the  most  significant  cause 
of  experimental  error,  that  due  to  uncertainty  in  the  separation 
measurements. 

The  sediment  is  disturbed  each  time  a  transducer  is 
repositioned  to  take  travel  time  and  amplitude  measurements  at  a  new 
separation.  It  is  also  disturbed  after  each  set  of  measurements,  when 
the  pore  fluid  is  circulated  and  glycerine  is  added  to  it  to  modify  its 
properties.  The  sample  is  also  vibrated  during  this  process,  so  it  is 
disturbed  even  more.  Any  such  disturbance  of  the  sediment  would  affect 
the  packing  of  the  indivioual  grains.  There  has  been  some  concern  that 
the  acoustic  propagation  through  sediments  was  very  dependent  upon  this 
packing.  If  the  propagation  were  very  dependent  on  the  packing,  it  would 
be  very  difficult  to  repeat  the  results.  The  repeatability  of  our  results 
suggests  that  the  propagation  is  not  very  dependent  on  the  packing  of  the 
individual  grains,  at  least  in  the  case  of  compressional  wave 
propagation.  This  is  fortunate,  because  it  would  be  very  difficult  and 
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tedious  to  control  the  pocking  of  gloss  beads  and  sediments  of  sizes  this 
small.  The  good  correlation  between  the  data  and  their  leest  squares 
fits.,  shown  in  Figs.  8  through  11,  also  supports  this  conclusion.  We 
conjecture  that  the  actual  pocking  of  the  glass  beads,  or  sediment,  would 
have  more  of  an  effect  on  shear  wave  results. 

Problems  associated  with  the  coupling  between  the  transducers 
and  the  sediment  could  also  affect  the  repeatability  of  the  results.  When 
there  is  good  coupling  between  the  transducer  and  the  sediment,  a  large 
amount  of  energy  is  transferred  from  the  transducer  to  the  sediment  in 
the  case  of  the  transmitter,  end  vice-versa  for  the  receiver.  Due  to  a 
difference  of  impedances  between  the  piezoelectric  element  and  the 
sediment,  all  of  the  energy  will  not  be  transmitted,  but  if  the  coupling  is 
good,  the  amount  of  energy  transferred  will  be  a  maximum.  If  there  is 
not  good  coupling  between  the  transducers  arid  the  sediment,  the 
measured  voltage  amplitude  will  be  less  then  the  amplitude  of  a  signal 
produced  by  gcod  coupling.  Since  the  radiating  surface  of  our 
transducers  was  ver  y  much  larger  than  the  diameter  of  the  glass  beads, 
the  coupling  betv/een  the  transducers  and  the  sediment  was  always  good 
and  the  energy  transferred  from  the  transducer  to  the  sediment,  and 
vice-versa,  was  always  adequate.  Several  observations  support  this 
conclusion.  First,  a  strong  signal  was  received  immediately  after  the 
transducer  was  inserted.  Nevertheless,  the  entire  sample  was  gently 
vibrated  for  twenty  to  thirty  minutes,  and  then  the  sediment  set 
undisturbed  for  another  twenty  to  thirty  minutes  more  before  a 


measurement  was  taken.  The  signal  was  monitored  for  at  least  three 
more  hours  before  the  transducer  was  repositioned.  During  this  period, 
the  time  delay  end  amplitude  readings  varied  only  by  amounts  attributed 
to  instrumental  uncertainties  The  good  correlation  between  the  data 
end  the  least  square  fits  also  supports  the  conclusion  that  there  was  no 
appreciable  variation  in  the  coupling  between  the  transducer  and  the 
sediment  during  the  period  the  readings  were  being  taken. 

In  the  discussion  of  experimental  errors  in  the  phase  velocity 
and  attenuation,  no  mention  has  been  made  of  the  accuracy  with  which 
the  viscosity  was  determined.  As  mentioned  in  Section  1 1  I  D.  on  the 
preparation  of  the  sediment,  the  viscosity  was  determined  by  measuring 
the  specific  gravity  of  the  pore  fluid  using  a  hydrometer,  end  then 
reading  the  value  of  viscosity  from  a  table  containing  the  specific 
gravity  end  viscosity  of  varying  concentrations  of  glycerine-water 
solutions.  The  error  in  the  viscosity  is  difficult  to  estimate,  host  of 
this  error  is  believed  to  be  due  to  the  variations  of  the  temperature  in 
the  room  where  the  experiment  was  performed.  The  temperature  in  the 
laboratory  was  usually  between  24°C  and  26°C,  so  most  pore  fluid 
viscosities  were  determined  assuming  that  the  room  temperature  was 
25°C.  The  theoretical  predictions  were  made  assuming  that  the 
temperature  was  25C'C.  For  the  most  part  this  produced  acceptable 
results,  especially  for  the  triangle  date  points  of  Fig.  15.  Luckily,  the 
laboratory  temperature  stayed  consistently  between  24°C  and  26°C 
during  the  time  these  data  points  were  taken  . 


For  the  other  measurements,  the  temperature  varied  between 
22°C  and  26°C.  This  4°  variation  in  temperature  can  produce  a  change  in 
viscosity  of  over  105?.  Not  only  did  this  affect  the  accuracy  at  which  the 
viscosity  was  determined,  but  it  also  affecied  the  velocity  and 
attenuation  results,  which  are  temperature  and  viscosity  dependent. 
Therefore  temperature  variation  could  be  a  source  of  some  of  the 
uncertainties  in  the  attenuation  and  velocity  measurements.  In  future 
investigations,  it  would  be  advisable  to  perform  similar  experiments  in 
constant  temperature  environments. 

Taking  the  experimental  uncertainties  that  have  been  discussed 
into  account,  the  B;ot-Stoll  theory  appears  to  adequately  describe  the 
effects  of  the  fluid  properties  on  acoustic  wave  propagation  in  saturated 
sediments,  at  least  for  compressionel  waves  of  the  first  type.  Similar 
experiments  investigating  shear  wave  propagation  and  the  compressionel 
wave  of  the  second  type  need  to  be  performed.  We  observe  that  these 
types  of  measurements  are  more  difficult  to  make. 


V!.  DISCUSSION  AND  CONCLUSIONS 


The  purpose  of  the  work  presented  here  wos  to  design  and 
perform  an  experiment  that  would  examine  the  dependence  of 
compressionel  wave  propagation  in  fluid  saturated  sands  on  the 
properties  of  the  pore  fluid.  Since  this  dependence  was  predicted  by  the 
Biot-Stoll  equations,  it  was  hoped  that  these  measurements  would  be  a 
critical  test  of  the  theory.  Experiments  of  this  type  had  already  been 
performed  by  Shirley  (Bedford,  Appendix  A,  1982),  and  then  later  by  Elliot 
(Bedford,  1982).  However,  due  to  large  scatter  in  the  experimental  date, 
no  definite  conclusions  could  be  made  about  their  results. 

The  measurement  techniques  used  in  this  investigation  were 
refinements  of  techniques  developed  and  used  by  Shirley  (Shirley,  1977; 
Bedford,  1982),  Bell  (Bell,  1979),  and  Elliot  (Bedford,  1982).  All  three  of 
these  investigators  measured  the  phase  velocity  and  attenuation  of  both 
compressionel  and  shear  waves.  After  making  a  few  modifications  in 
their  technique,  new  measurements  have  been  taken  and  shown  to  be 
repeatable.  The  main  improvement  over  the  earlier  work  resulted  from 
redesigning  the  compressionel  wave  transducers  so  that  the  geometrical 
spreading  end  directivity  could  be  predicted.  With  th*s  improvement, 
spherical  spreading  was  established  at  the  separations  at  which  these 
measurements  were  taken.  Thus,  the  losses  due  to  geometrical  spreading 
were  accurately  accounted  for  when  analyzing  the  data.  It  was  also 
possible  to  make  reasonable  estimates  of  the  experimental 
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uncertainties.  Because  of  the  paramount  difficulties  involved  in  making 
repeatable  shear  wave  measurements,  these  measurements  were 
postponed. 

Taking  the  experimental  uncertainty  into  account,  the  data  are 
in  good  agreement  with  the  theory.  Since  parameters  were  chosen  so 
that  the  theoretical  values  were  in  agreement  with  the  experimental 
results  for  the  case  of  pure  water,  the  present  results  indicate  that  the 
theory  predicts  the  same  trend  in  attenuation  and  phase  velocity,  with 
changing  properties  of  the  pore  fluid,  that  is  exhibited  by  the  data.  But 
since  there  was  also  good  agreement  with  the  fresh  water  velocity 
results  of  Shirley,  taken  at  approximately  the  same  depth,  the  choices  of 
the  real  parts  of  the  bulk  and  shear  moduli  appear  to  be  reasonable 
estimates  of  these  parameters.  Further  experimental  investigations 
examining  the  depth  and  frequency  dependence  of  compressional  wave 
propagation  in  sedimenis  ere  needed  to  support  this  conclusion. 

These  compressional  wave  results,  on  the  other  hand,  do  not 
give  a  strong  indication  of  the  values  of  the  imaginary  parts  of  these 
moduli.  To  get  more  reliable  estimates  of  these  parameters,  shear  wave 
measurements  in  the  seme  size  glass  beads  need  to  be  token.  Shear  wove 
velocity  measurements  ore  easily  accessible  with  only  a  modest 
investment  in  time  and  equipment. 

On  the  other  hand,  it  is  believed  that  shear  wav?  attenuation 
measurements  would  be  much  more  difficult.  One  of  these  difficulties, 
especially  in  the  gloss  beads,  would  be  the  coupling  between  the 
transducer  and  the  sediment.  Once  this  problem  has  been  surmounted, 


establishing  the  directivity  of  shear  wave  transducers  would  also  be 
tedious  and  time  consuming.  However,  if  a  procedure  similar  to  the  one 
described  here  were  used,  this  would  be  an  essential  step.  Only  in  this 
way  could  an  experimenter  be  sure  that  the  amplitude  decay  being 
measured  was  due  to  the  increased  separation  of  the  transducers,  and  not 
their  directivity. 

It  would  also  be  important  to  establish  which  pert  of  the  total 
losses  would  be  due  to  the  geometrical  spreading  of  the  transducers. 
This  would  probably  be  the  most  difficult  of  the  three  tests  to  conduct. 
The  most  common  type  of  transducers  used  in  experimental  shear  wave 
propagation  studies  are  shear  wave  bender  elements.  The  problem  would 
be  to  locate  a  lossless,  elastic  medium  in  which  the  bender  elements 
would  work  so  tnat  the  losses  due  to  the  spreading  of  the  signal  could  be 
determined. 

Fortunately  for  our  case,  Rayleigh  had  already  solved  for  the 
ferfield  pressure  produced  by  a  moving  piston  in  a  rigid  baffle.  All  that 
we  needed  to  do  was  verify  that  this  model  adequately  predicted  the 
measured  characteristics  of  our  transducers  end  that  our  measurements 
were  taken  in  the  farfield.  It  would  be  helpful  if  a  similar  model  could 
be  derived  to  describe  the  stress  field  or  the  displacement  field  of  shear 
bender  elements.  With  a  model  of  this  type,  establishing  the 
characteristics  of  this  type  of  transducer  would  be  simpler. 
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APPENDIX  A 


COMPUTER  PROGRAM 


The  FORTRAN  program  presented  here  is  a  simple  algorithm  that 
calculates  the  phase  velocity  and  attenuation  of  sound  in  saturated 
sediments  using  the  Biot-Stoll  theory.  The  program  solves  the  quadratic 
formula  (11-3),  from  Chapter  II  for  different  values  of  viscosity,  density, 
ond  bulk  modulus  of  the  pore  fluid.  The  fluid  properties  used  in  this 
study  are  shown  in  Table  VI  of  Chapter  V  and  were  fed  to  the  program  by 
the  input  data  file  DATAIN.  The  physical  properties  of  the  solid  and  solid 
matrix  were  either  entered  interactively  or  in  the  DATA  statement.  The 
values  of  these  properties  are  listed  in  Table  V  of  Chapter  V. 

The  experimental  results  of  Table  VII  were  also  entered  into  the 
program  through  DATAIN.  The  program  plots  the  theoretical  arid 
experimental  results,  shown  in  Fig.  15,  so  that  they  can  be  compared.  A 
list  of  the  parameters  used  in  the  program,  along  with  their 
corresponding  FORTRAN  and  mathematical  symbols,  is  given  in  Table  VIII, 
followed  by  a  listing  of  the  FORTRAN  program  BIVO. 
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Table  VIII:  Parameters  Used  in  FORTRAN  program  BIVO. 


PARAMETER  MATH  SYMBOL 

Pore  Size  ap 


Virtual  Drag  Coefficient 

Virtual  Mass  Coefficient 

Biot  Coefficients 
(complex) 

(real) 

Permeability 
Grain  Diameter 

Porosity 
Solid  Density 

Fluid  Density 

Mixture  Density 

Viscosity 


b 

c 

C 

H 

M 

D 

Bo 

'Vi 

♦ 


Ps 

pf 

p 


Pc 


'n 


Phase  Velocity 


co 


FORTRAN  SYMBOL 
A 

B 

C 

CB 

HB 

MB 

D 

BO 

DM 

PHI 

ROS 


ROF 


RO 


ROC 

ETA,  VISC 
VISC2,  VI3C3 


CPHASE 
Cl,  C2,  C3 
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PARAMETER 


MATH  SYMBOL 


FORTRAN  SYMBOL 


Attenuation 

OC 

ATTEN,  ALPHA 
ALPHA2,  ALPHA 

Frequency 

f 

FR 

Angular  Frequency 

<0 

OMEGA 

Fluid  Bulk  Modulus 

Kf 

KBF 

Solid  Bulk  Modulus 

Ks 

KBS 

Bulk  Modulus  of  Brained 

Solid  Matrix  (complex) 

KBB 

Shear  Modulus  of  Drained 

Solid  Matrix  (complex) 

MUB 

Compressional  Logarithmic 
Decrement 

£c 

DELC 

Shear  Logarithmic  Decrement 

6s 

DELS 
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APPENDIX  B 


EXPERIMENTAL  DATA 


BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


AMPLITUDE  -  volts 

TIME  -  MS 

SEPARATION 

0.3720 

65.88 

5.65 

0.2690 

77.63 

7.65 

0.2004 

88.98 

9.65 

0.1570 

100.12 

11.65 

0.1300 

112.07 

13.65 

0.1074 

121.38 

15.65 

0.0892 

133.23 

17.65 

0.0777 

142.31 

19.65 

VELOCITY:  1020.0  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  23.6  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9995 

ATTENUATION:  2.38  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.4?  Np/m 


BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


AMPLITUDE  -  volts 

TIME  -  jis 

SEPARATION  -  cm 

0.1286 

62.21 

5.65 

0.0907 

73.60 

7.65 

0.0690 

84.04 

9.65 

0.0542 

94.83 

1 1.65 

0.0440 

105.90 

13.65 

0.0366 

1 15.97 

15.65 

0.0309 

126.98 

17.65 

0.0266 

137.16 

19.65 

VELOCITY:  1869.7  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  9.0  rn/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 


ATTENUATION:  2.36  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.12  Np/m 
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BEAD  TYPE:  MS- 

-MH 

FREQUENCY: 

105.0  kHz 

VISCOSITY:  1.07  cP 

AMPLITUDE  -  volts 

TIME  -  MS 

SEPARATION  -  cm 

0.1310 

61.69 

5.65 

0.0912 

72.61 

7.65 

0.0696 

82.62 

9.65 

0.0547 

93.72 

1 1.65 

0.0441 

104.37 

13.65 

0.0372 

115.05 

15.65 

0.0310 

126.53 

17.65 

0.0268 

136.68 

19.65 

VELOCITY:  1861.3  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  6.8  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 

ATTENUATION:  2.46  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.25  Np/m 
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BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  1.22  cP 


AMPLITUDE  -  volts 

TIME  -  ps 

SEPARATION 

0.1325 

61.17 

5.65 

0.0928 

71.66 

7.65 

0.0702 

82.34 

9.65 

0.0550 

92.60 

1  1.65 

0.0448 

103.58 

13.65 

0.0370 

1 14.31 

15.65 

0.0321 

124.96 

17.65 

0.0270 

135.59 

19.65 

VELOCITY:  1878.9  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  2.6  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 

ATTENUATION:  2.42  Np/m 

5TANDARD  DEVIATION  OF  ATTENUATION:  0.36  Np/m 
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BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  1.33  CP 


AMPLITUDE  -  volts 

TIME  -  jiS 

SEPARATION 

0.1280 

60.74 

5.66 

0.0916 

71.00 

7.66 

0.0697 

81.30 

9.63 

0.0544 

92.02 

11.65 

0.0445 

102.42 

13.62 

0.0375 

1 12.98 

15.61 

0.0320 

123.17 

17.61 

0.0274 

134.40 

19.63 

VELOCITY:  1898.6  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  5.7  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 

ATTENUATION:  2.18  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.22  Np/m 


BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  1.33  cP 


AfIPLitUDfc  -  volts 

tint  -  jjs 

SEPARATION 

0.1310 

60.97 

5.67 

0.0920 

71.20 

7.64 

0.070? 

81.59 

9.62 

0.0546 

92.60 

l  1.64 

0.0443 

103.14 

13.63 

0.0370 

1 13.4? 

15.62 

0.0316 

124.25 

17.65 

0.0275 

134.22 

19.58 

VELOCITY:  1893.3  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  4.2  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 

ATTENUATION:  2.35  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.33  Np/m 


BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  0.893  cP 


AMPLITUDE  -  volts 

TIME  -  us 

SEPARATION  -  cm 

0.1193 

63.61 

6.0 

0.0890 

73.32 

8.0 

0.0654 

85.83 

10.0 

0.0526 

95.53 

12.0 

0.0430 

106.16 

14.0 

0.0352 

117.32 

16.0 

0.0300 

128.65 

18.0 

0.0258 

139.50 

20.0 

VELOCITY:  1839.4  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  15.2  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9998 

ATTENUATION:  2.50  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.29  Np/m 


BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  1.15  cP 


AMPLITUDE  -  volts 

TIME  -  MS 

SEPARATION  -  cm 

0.1218 

62.70 

6.1 

0.0866 

73.42 

8.1 

0.0654 

84.32 

10.1 

0.0516 

94.61 

12.0 

0.042? 

104.86 

14.0 

0.0354 

1 15.86 

16.0 

0.0299 

126.78 

16.0 

0.0256 

136.53 

20.0 

VELOCITY:  1873.8  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  b.7  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 

ATTENUATION:  2.66  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.25  Np/m 


BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kH2  VISCOSITY:  1.48  cP 

AMPLITUDE  -  volts  TIME  -  jis  SEPARATION  -  crn 


0.1  190 

62.02 

6.05 

0.0850 

72.65 

8.09 

0.0646 

82.73 

10.00 

0.0518 

92.36 

12.05 

0.0416 

102.45 

14.02 

0.0340 

112.99 

16.03 

0.0280 

123.55 

18.05 

0.0238 

132.91 

20.00 

VELOCITY:  1966.7  m/s 

STANDARD 

DEVIATION  OF  VELOCITY: 

9.4  m/s 

LINEAR  CORRELATION  COEFFICIENT:  0.9999 
ATTENUATION:  2.99  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.41  Np/m 
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BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  1.74  cP 

AMPLITUDE  -  volts  TIME  -  pis  SEPARATION  -  cm 


0.1  193 

61.24 

6.06 

0.0870 

71.06 

8.02 

0.0649 

81.55 

10.02 

0.0514 

91.70 

12.05 

0.0420 

102.19 

14.05 

0.0350 

112.25 

16.02 

0.0304 

121.78 

18.02 

0.0258 

131.66 

20.00 

VELOCITY:  1973.3  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  8.9  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 

ATTENUATION:  2.41  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.33  Np/m 
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BEAD  TYPE;  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  2.25  cP 


AMPLITUDE  -  volts 

TIME  -  jis 

SEPARATION  -  cm 

0.1180 

61.42 

6.06 

0.0860 

70.69 

8.03 

0.0645 

80.80 

10.06 

0.0515 

90.20 

12.00 

0.0419 

100.52 

14.04 

0.0348 

110.25 

16.03 

0.0298 

119  87 

18.00 

0.0258 

129.16 

20.00 

VELOCITY:  2044.1  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  7.3  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 

ATTENUATION:  2.42  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.19  Np/m 
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BEAD  TVPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  2.59  cP 


AMPLITUDE  -  volts 

TIME  -  MS 

SEPARATION  -  cm 

0.1180 

60.55 

6.02 

0.0843 

70.45 

8.03 

0.0649 

79.30 

10.01 

0.0500 

89.72 

12.04 

0.0410 

98.97 

13.99 

0.0336 

109.21 

16.04 

0.0285 

118.76 

18.03 

0.0250 

127.41 

19.85 

VELOCITY:  2067.5  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  9.2  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 


ATTENUATION:  2.67  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.24  Np/m 


BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  2.81  cP 


AMPLITUDE  -  volts 

TIME  -  JJS 

SEPARATION  -  cm 

0.1 190 

60.41 

6.04 

0.0846 

70.00 

8.03 

0.0640 

79.91 

10.05 

0.0505 

89.64 

12.02 

0.0408 

99.09 

14.01 

0.0340 

108.63 

16.00 

0.0282 

119.06 

18.06 

0.0245 

126.97 

19.96 

VELOCITY:  2074.5  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  13.0  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 

ATTENUATION:  2.80  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.16  Np/m 
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BEAD  TYPE:  MS-MH 

FREQUENCY:  105.0  kHz  VISCOSITY:  3.23  CP 


AMPLITUDE  -  volts 

TIME  -  jjs 

SEPARATION  -  cm 

0.1160 

60.36 

6.07 

0.0830 

69.85 

8.06 

0.0640 

79.20 

10.06 

0.0505 

80.30 

12.01 

0.0410 

98.12 

14.03 

0.0336 

107.62 

16.02 

0.0280 

1 17.86 

18.01 

0.0240 

126.48 

20.05 

VELOCITY:  2089.8  m/s 

STANDARD  DEVIATION  OF  VELOCITY:  12.1  m/s 
LINEAR  CORRELATION  COEFFICIENT:  0.9999 


ATTENUATION:  2.80  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.26  Np/m 
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BEAD  TYPE:  MS-MH 


FREQUENCY: 

105.0  kH2 

VISCOSITY:  2.46  CP 

AMPLITUDE  -  volts 

TIME  -  ps 

SEPARATION  -  cm 

0.1 158 

61.06 

6.07 

0.0820 

71.26 

8.09 

0.0616 

80.76 

10.07 

0.0484 

90.82 

1 2.06 

0.0.396 

100.46 

14.03 

0.0328 

1  10.37 

16.02 

0.0274 

1  19.43 

18.01 

0.0231 

129.68 

19.98 

VELOCITY:  2034.7  m/s 
STANDARD  DEVIATION  OF  VELOCITY:  8.2  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 

ATTENUATION:  2.99  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.19  Np/m 


BEAD  TYPE: 

FREQUENCY:  105.0  kH2 


MS-MH 

VISCOSITY:  2.09  cP 


AMPLITUDE  -  volts 


TIME  -  jis  SEPARATION  -  cm 


0.1210 

0.0870 

0.0660 

0.0526 

0.0429 

0.0355 

0.0302 

0.0262 


61.73 

6.01 

71.80 

8.02 

82.32 

10.06 

91.83 

12.02 

101.86 

14.00 

1 12.19 

16.02 

122.23 

18.02 

132.08 

19.99 

VELOCITY:  1986.6  m/s 
STANDARD  DEVIATION  OF  VELOCITY;  3.6  m/s 
LINEAR  CORRELATION  COEFFICIENT:  1.0000 


ATTENUATION:  2.39  Np/m 

STANDARD  DEVIATION  OF  ATTENUATION:  0.19  Np/m 


.  v j v r ”»  w .  ’  -  ’ 


$ 

Ti 


REFERENCES 


Bedford,  A.,  G.  Ellis,  K.  P.  Elliot,  S.  K.  Mitchell,  D.  J.  Shirley, 
"Acoustical  Properties  of  Sediments,"  Applied  Research 
Laboratories  Technical  Report  No.  82-47  (ARL-TR-82-47),  Applied 
Research  Laboratories,  The  University  of  Texas  at  Austin  (1982). 

Bedford,  A.,  R.  D.  Costley,  and  M.  Stern,  "On  the  Drag  and  Virtual  Mass 
Coefficients  in  Biot’s  Equations,”  J.  Acoust.  Soc.  Am.  76,  1804-1809 
(1984). 

Bell,  David  W.,  "Shear  Wave  Propagation  in  Unconsolidated  Fluid 
Saturated  Porous  Media,”  Applied  Research  Laboratories  Technical 
Report  No.  79-31  (ARL-TR-79-31),  Applied  Research  Laboratories, 
The  University  of  Texas  at  Austin  (1979). 

Bevington,  Philip  R.,  “Data  Reduction  and  Error  Arialusis  for  the 
Physical  Sciences”  (McGraw-Hill,  Book  Co.,  Inc.,  New  York,  1969). 

Biot,  M  A.,  “Theory  of  Propagation  of  Elastic  Waves  in  a  Fluid-Saturated 
Porous  Solid.  I.  Low  Frequency  Range,”  J.  Acoust.  Soc.  Am.  28. 
168-178  (1956a). 

Biot,  M.  A.,  "Theory  of  Propagation  of  Elastic  Waves  in  a  Fluid-Saturated 
Porous  Solid.  II.  Higher  Frequency  Range,"  J.  Acoust.  Soc.  Am  28, 
179-191  (1956b). 

Biot,  M.  A.,  ”  Generalized  Theory  of  Acoustic  Propagation  in  Porous 
Dissipative  Media,"  J.  Acoust.  Soc.  Am.  34,  1 254- 1 264  ( 1 956c). 

Bleckstock,  David  T.,  Notes  from  "Acoustics  I  and  II,"  courses  given 
in  the  Departments  of  Mechanical  Engineering  and  Electrical 
Engineering  at  The  University  of  Texas  at  Austin. 


95 


Greenspan.,  Martin,  "Acoustic  Properties  of  Liquids. “  in  American 
Institute  of  Phusics  Handbook.  3rd  ed.,  Dwight  E.  Gray, 
ed.  (McGraw-Hill  Book  Co.,  Inc.,  New  York,  1963),  pp.  3-79. 

Hovem,  J.  M.,  and  G.  D.  Ingram,  "Viscous  Attenuation  of  Sound  in 
Saturated  Sand,"  J.  Acoust.  Soc.  Am.  66,  1807-181 1  (1979). 

Kinsler,  L.  E.,  A.  R.  Frey,  A.  B.  Coppens,  and  J.  V.  Sanders,  Fundamentals 
of  Acoustics.  3rd  ed.  (John  Wiley  &  Sons,  Inc.,  New  York,  1982), 
p.  461. 

Sheely,  Madison  L.,  Industrial  and  Engineering  Chemistry  24, 
1060-64(1932). 

Shirley,  Donald  J.,  and  David  W.Bell,  "Acoustics  of  in  Situ  and 
Laboratory  Sediments,"  Applied  Research  Laboratories  Technical 
Report  No.  78-36  (ARL-TR-78-36),  Applied  Reseorch  Laboratories, 
The  University  of  Texas  at  Austin  (1978). 

Stoll,  R.  D.,  "Acoustic  Woves  in  Soturoted  Sediments,"  in  Phusics  of 
Sound  in  Marine  Sediments.  L.  D.  Homptori  ,  ed.  (Plenum  Press, 
New  York,  1974),  pp.  19-39. 

Stoll,  R.  D.,  "Acoustic  Woves  in  Oceon  Sediments,"  Geophysics  42, 
715-725  (1977). 


WCTCTY?. ers’I  i’.,'»S7"*"1-v'irT.,V'  v  -'n.-vr--,^' *:■*  ,T*  tjv 


ACKNOWLEDGMENTS 

The  author  was  fortunate  to  receive  contributions  and  support 
from  several  persons.  Dr.  A.  Bedford  supervised  this  research  and  was  an 
appreciated  source  of  encouragement  and  patience.  Dr.  Paul  Vidmar  made 
many  contributions  and  suggestions. 

Dr.  Ilene  Busch-Vishniac  answered  my  many  questions, 
particularly  those  involving  the  design  of  the  transducers  and  error 
analysis.  The  transducer  design  was  primarily  the  result  of  two 
conversations  with  Dr.  Elmer  Hixson.  Mr.  Ken  Kirksey  made  several 
contributions  in  the  development  of  the  transducers,  particularly  in  their 
construction.  Dr.  David  Blackstock  was  very  encouraging  and  referred  me 
to  several  helpful  sources.  I  was  also  fortunate  to  visit  Mr.  James 
Matthews'  lab  at  NCiRDA  and  to  witness  some  of  the  progress  he  has  made 
in  experimental  sediment  acoustics.  Dr.  George  B.  Thurston  advised  me  on 
the  viscosity  measurements.  I  am  indebted  to  Mr.  Donald  J.  Shirley  and 
Mr.  Kenneth  P.  Elliot,  who  both  preceded  me  in  the  sediments  lob  at 
ARL:UT.  The  experiment  was  the  brain  child  of  Mr.  Shirley. 

Ms.  Jo  Liridberg  and  Mr.  James  TenCate  made  several 
contributions  concerning  instrumentation. 

This  research  was  supported  by  the  Office  of  Navel  Research, 
Code  425GG  and  Code  400,  end  was  submitted  to  the  Graduate  School  of 
The  University  of  Texas  at  Austin  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Master  of  Science  in  Engineering. 


97 


.'•'.'•■/.•.'■.•.'•'.V.V.V-.V.V.'--.  O'.Vo.'-V 


\r*  ,-T*  *-#>^1  ar.  '*j!t«  *^T£z?.zr 


8 


s 

.** 


Copy  No, 


2 

3 

4 

5 

6 

7 

8 


9 

10 

11 


12 

13 

14 


15 


16 


May  1985 


DISTRIBUTION  LIST  FOR 
ARL-TR-85-14 

UNDER  CONTRACT  N00014-80-C-0490 


Office  of  Naval  Research 
495  Summer  Street 
Boston,  MA  02210 
Attn:  R.  L.  Sternberg 


Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  MS  39529 
Attn:  E.  D.  Chaika  (Code  530) 

D.  B„  King  (Code  321) 

R.  Wagstaff  (Code  245) 

K.  Gilbert  (Code  221) 

A  Upn  7p1 

G.*  Morris  (Code  425GG) 

J.  Matthews  (Code  362) 


Chief  of  Naval  Research 
Department  of  the  Navy 
Arlington,  V A  22217 

Attn:  M.  McKisic  (Code  425-OA) 

R.  Obrochta  (Code  425-AR) 

F.  E.  Saal field  (Code  400) 


Director 

Naval  Research  Laboratory 
Department  of  the  Navy 
Washington,  DC  20375 
Attn:  0.  Diachok  (Code  5128) 
R.  Dicus  (Code  5128) 

R.  Gragg 


Lawrence  Livermore  National  Laboratory 
P.  0.  Box  808,  L-200 
Livermore,  CA  94550 
Attn:  J.  G.  Berryman 


Planning  Systems,  Inc. 
7900  Westpark  Drive 
Suite  600 
McLean,  VA  22102 
Attn:  B.  Brunson 


99 


PV 


Distribution  List  for  ARL-TR-85-14  under  Contract  N00014-80-C-04S0 
(cont'd) 

Copy  No. 

Lamont-Daherty  Geological  Observatory 
Columbia  Uriversit.y 
Palisades,  NY  10964 

17  Attn:  R.  D.  Stoll 

18  A.  Mark  Bedford,  ARLtUT 

19  H.  Boehme,  ARL :UT 

,!0  Ilene  J,  Busch-Vishniac ,  ARL:UT 

21  Karl  C.  Focke,  ARL : UT 

22  Thomas  A.  Griffy,  ARLtUT 

23  John  M.  Huckabay,  ARL:UT 

24  Robert  A.  Koch,  ARL :UT 

25  Stephen  K.  Mitchell,  ARLtUT 

26  T.  G.  Muir,  ARLtUT 

27  David  W.  Oakley,  ARLtUT 

28  Clark  S.  Penrod,  ARLtUT 

29  Morris  Stern.  ARLtUT 

30  James  A.  TenCate,  ARLtUT 

31  Paul  J.  Vidmar,  ARLtUT 

32  Gary  R.  Wilson,  ARLtUT 

33  H.  Yew,  ARLtUT 
Library,  ARLtUT 


34 


